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ABSTRACT 


A thermohydrodynamic analysis is presented and a computer code developed 
for prediction of the static and dynamic force response of hydrostatic journal 
bearings (HJBs), annular seals or damper bearing seals, and fixed arc pad bearings 
for cryogenic liquid applications. The study includes the most important flow 
characteristics found in cryogenic fluid film bearings such as flow turbulence, fluid 
inertia, liquid compressibility and thermal effects. The analysis and computational 
model devised allow the determination of the flow field in cryogenic fluid film 

bearings along with the dynamic force coefficients for rotor-bearing stability 
analysis. 


Journal model: ASME Journal of Tribology. 
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EXECUTIVE SUMMARY 

A thermohydrodynamic analysis and computer code for the prediction of the 
static and dynamic force response of fluid film bearings for cryogenic applications 
have been completed on 1993. The current technological needs require reliable 
and resilient fluid film bearing designs to provide maximum operating life with 
optimum rotordynamic characteristics at the lowest cost. The analysis and code 
constitute practical tools for the prediction of performance and design of cryogenic 
liquid hydrostatic journal bearings, annular pressure seals or damper bearings, and 
cylindrical pad hydrodynamic bearings. 

The motion of a cryogenic liquid on the thin film annular region of a 
fluid film bearing is described by a set of mass, momentum conservation, and 
energy transport equations for the primitive turbulent bulk-flow variables, and 
accompanied by thermophysical state equations for evaluation of the fluid material 
properties. Zeroth-order equations describe the fluid flow field for a journal static 
equilibrium position, while first-order linear equations govern the fluid flow for 
small amplitude journal center translational motions and journal axis conical 
motions. Solution to the zeroth-order flow field equations provides the bearing 
flow rate, load capacity, restoring moments and torque. Solution to the first- 
order equations determines the rotordynamic force and moment coefficients due 
to journal lateral and angular motions. The analysis includes the effects of flow 
turbulence, fluid inertia, liquid compressibility and thermal energy transport on 
the performance of cryogenic liquid bearings. 

The numerical predictions from the computational program ( hydrosealt ) 
correlate favorably with experimental results available in the literature. Numerous 
examples are provided to demonstrate the effectiveness of the program to analyze 
cryogenic fluid film bearings as well as conventional viscous lubricant bearings. 
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INTRODUCTION 

Hydrostatic Journal Bearings (HJBs) are the ideal candidates to replace roller 
bearings as support elements in cryogenic turbomachinery. These bearings along 
with hydrostatic annular seals will be used for primary space-power applications 
due to their long lifetime, low friction and wear, significant load capacity, and 
large direct stiffness and damping force coefficients. HJBs, unlike rolling element 
bearings, have no DN limit. Rotating machinery free of this constraint can 
operate at larger speeds with better efficiency and reduced overall weight and 
size. Durability in HJBs is assured by the absence of contact between stationary 
and moving parts during steady-state operation, while long life reduces the 
frequency of required overhauls. Despite these attractive features, fluid film 
bearing stability considerations and thermal phenomena along with two-phase 
flow operation are a primary concern for high speed operation with large pressure 
differentials. Fluid film bearing stability is essentially related to hydrodynamic 
and liquid compressibility effects, while thermal effects are of importance due to 
flow turbulence with substantial energy dissipation. 

A comprehensive literature review on the subject of hydrostatic bearings and 
annular seals relevant to cryogenic liquid applications has been presented elsewhere 
(San Andres, 1990a-c, 1991a-c, 1992a, 1993a). A brief discussion of the relevant 
literature pertaining to the analysis, design and testing of fluid film bearings for 
cryogenic applications over the last few years follows. 

Reddecliff and Vohr (1969) initially studied HJBs for use in high-pressure 
cryogenic rocket engine turbopumps. The inertia effect at the edges of the recesses 
changed the pressure distribution, reduced the flow rate but did not affect the 
total bearing load capacity. The non linear fluid advective inertial terms could 
not be accounted for due to numerical difficulties. Variable fluid properties were 
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treated as linear between those at the supply and discharge pressures, and steady- 
state predictions were reported to agree well with experimental results. The 
scatter in the measured recess pressures was attributed to variations in bearing 
clearance caused by temperature differences between the bearing and shaft and 
to measurement inaccuracies. 

Artiles et al. (1982) presented a numerical solution to the static and dynamic 
performance characteristics of hydrostatic journal bearings. A turbulent Reynolds 
equation with constant fluid properties was solved by the column-matrix method, 
while a Newton-Raphson scheme was implemented for efficient calculation of the 
recess pressures. Turbulent-to-laminar flow power-loss ratios were reported to 
be in the range of 25 to 30 for the Reynolds numbers considered. Even though 
there were neither energy considerations nor thermal effects in the analysis, large 
temperature rises (up to 24.5°(7) in the fluid film were reported for LO 2 bearings, 
while temperature rises in LH 2 bearings were negligibly small. 

Braun et al.( 1987a, b) introduced a comprehensive THD analysis for cryogenic 
liquid hydrodynamic and hydrostatic journal bearings. On the fluid film region, 
a variable-properties Reynolds equation was coupled to a 2-D energy transport 
equation. The heat transfer to the bounding solids (shaft and bush) was analyzed 
in its three dimensional complexity. Bulk-flow heat transfer coefficients were used 
to represent the boundary conditions at the fluid/solid interfaces, and fluid inertia 
effects were considered only at the pocket’s edges with no recess volume-liquid 
compressibility effects. Braun et al.’s analysis regard the fluid flow as laminar, 
although large pressure differentials and rotational speeds were considered in the 
applications studied. For a liquid hydrogen HJB, a small temperature increment 
was found in the fluid film, and thermal effects were shown to be minimal 
relative to a constant properties liquid model. Calculations for a liquid oxygen 
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hydrodynamic journal bearing showed a large temperature rise (13 °K) for eccentric 
operation at 14,000 rpm with a profound impact on the mass flow rate and load 
capacity. The numerical predictions presented show circumferential flow Reynolds 
numbers as large as 100,000 with a laminar flow model. No definite conclusions 
in regard to the effect of heat transfer from the fluid film to the bounding solids 
are extracted from the analysis. 

San Andres (1990a) introduced the first full inertial, fully developed turbulent 
bulk-flow model for the analysis of incompressible liquids hydrostatic bearings. 
This model revealed the importance of flow turbulence and fluid inertia at the 
film lands and at the recess boundaries of typical high speed hydrostatic bearings. 
Results of the analysis show that a fluid inertia flow model when compared to a 
classical viscous model determines reduced leakage rates and hydrostatic forces 
while increasing the hydrodynamic force due to journal rotation. These effects 
then lead to the generation of larger direct damping and cross-coupled stiffness 
coefficients whith a reduction in the direct stiffness coefficients and the appearance 
of significant inertia force coefficients for liquids of large density such as water 
and liquid oxygen. San Andres (1991b) shows that moderate to large journal 
eccentricities have a pronounced effect on the force coefficients of HJBs with large 
hydrodynamic effects (high rotational speeds). For large journal center static 
excursions, the hydrodynamic effect produces force coefficients similar in form to 
those of a conventional circular journal bearing. Furthermore, orifice back-flow 
accompanied by a sudden drop on direct stiffness is likely to occur at large journal 
eccentricity operation. 

A study of the laminar flow field on a recess-land region of a typical 
high speed hydrostatic bearing geometry is given by San Andres and Velthuis 
(1992). Computational experiments revealed a substantial hydrodynamic pressure 
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generation on the downstream side of a recess followed by a sharp pressure drop 
at the recess edge next to a film region of small clearance. Calculated recess-edge 
entrance factors and recess shear coefficients form the basis of a simplified one- 
dimensional recess bulk-flow model implemented for the analysis of finite length 
ydrostatic bearings. Other recent analysis and numerical predictions studying 
the same flow geometry are given by Braun et al. (1993), and Hill et al. (1993). 

The Texas A&M hydrostatic bearing test facility directed has provided with 
a wealth of experimental data for the static and dynamic force characteristics of 
water lubricated hydrostatic bearings. Experimental measurements are routinely 
performed for bearings of different geometries and at journal speeds ranging from 
10,000 to 25,000 rpm and pressure supplies from 4 to 7 MPa. The test facility 
accomodates state of the art instrumentation and control, and possesses the most 
efficient and accurate bearing parameter identification method based on a real- 
time process and frequency domain algorithms. Kurtin et al. (1993) reported 
experimental data for the static performance characteristics of a 5 recess, water 
HJB for ,he °P erat ing conditions noted and three different bearing clearances. The 
experimental results correlate very well with predictions from the numerical model 
developed by San Andres (1990a, 1992a). It is noted also that accurate theoretical 
results depend greatly on the knowledge of the bearing operating clearance, and 
most importantly, on the orifice discharge coefficients. 

The stability characteristics of a simple rotor-bearing system are defined by 
the threshold speed of instability and the whirl frequency ratio. This instability is 
of the ’’hydrodynamic” type and solely due to the effect of journal rotational speed 
on the flow field of a bearing under dynamic journal motions. The threshold speed 
represents the journal rotational speed at which the bearings loose their effective 
damping and any small perturbation from an equilibrium position determines 
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unbounded rotor motions according to linear vibrations theory. The rotor-bearing 
first natural frequency (critical speed) depends on the bearing equivalent stiffness 
and the shaft stiffness. It is also well known that flexible rotors do present natural 
frequencies much lower than rigid rotors. The whirl frequency ratio denotes the 
ratio between the onset whirl frequency (typically the system first critical speed) 
and the threshold speed of instability. This ratio is independent of the flexibility 
of the rotating shaft and, it provides information on the stability characteristics of 
the fluid film bearings within the desired range of operation of a typical machine. 
Plain cylindrical journal bearings show a characteristic whirl frequency ratio equal 
to 0.50 for operation at small to moderate eccentricities (light loads) and hence 
denote that a simple rotor-bearing system will become unstable at a rotational 
speed equal to twice the system first critical speed. 

The results presented by San Andres (1990c, 1991b) have shown that, con- 
trary to generalized intuition, incompressible liquid hydrostatic bearings present 
a whirl frequency ratio identical to that of plain journal bearings. This condition 
then limits severely the application of HJBs to high speed, light weight turbo- 
machinery. Furthermore, cross-coupled inertia force coefficients for large density 
fluids have an adverse effect on the threshold speed of instability The experimen- 
tal results from the Hydrostatic Bearing Test Program at Texas A&M University 
verify closely the theoretical predictions and show, in some circumstances, the 
whirl frequency ratio to increase above 0.50 for low rotational speeds and large 
supply pressures (Franchek et ah, 1993, Mosher, 1993). 

Approximate theoretical models for prediction of the dynamic force response 
of hydrostatic bearings at the centered position have been reported also by San 
Andres (1990b, 1991a, 1992b). Approximate analytical solutions to the turbulent 
flow field bring insight into the understanding of the flow mechanics and provide 
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a fast and efficient tool for the preliminary design of hydrostatic journal bearings. 
Results of the analyses show the importance of recess volume fluid compress- 
ibility effects on the rotordynamic force response of typical HJB for cryogenic 
applications. In brief, HJBs handling highly compressible liquids such as LH 2 
for example, are prone to show a self-excited type instability of the ’’pneumatic 
hammer” type and could produce negative damping force coefficients for low fre- 
quency excitations. Dynamic operation under these conditions will then result 
in a whirl frequency ratio greater than 0.50, and consequently, it aggravates the 
’’hydrodynamic” instability problem since it reduces even further the safe range of 
operating speeds of a rotor-bearing system. This important result, although first 
reported by Redecliff and Vohr (1969), has been largely overlooked until recently. 

Hydrostatic bearings for cryogenic liquids of low viscosity require large levels 
of external pressurization to provide adequate load capacity and stiffness support. 
Typical pressure drops across a HJB can be as large as 30MPa and under 
these conditions the material properties of cryogenic liquids present substantial 
variations. These considerations lead San Andres (1992a) to develop a complex 
analysis for variable properties liquid HJBs. The fluid is regarded as barotropic 
with material properties, density and viscosity, as functions of a characteristic 
temperature and the local absolute pressure within the bearing flow region. Note 
that the designation of barotropic applies only to the fluid and not to the flow 
process within the fluid film bearing. The numerical model extracts the cryogenic 
liquid properties from accurate 32-term thermophysical state equations as given 
in the data base of McCarty(1986). Predictions of the isothermal analysis for 
a typical liquid hydrogen HJB show that the fluid properties affect considerably 
the static and dynamic force characteristics of the bearing and determine a whirl 
frequency ratio much larger than 0.50 for centered journal operation. This result 
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is a direct consequence the fluid compressibility effect at the recess volumes and at 
film lands. Sump pressures are relatively low at the discharge side of the bearing 
and determine small values of LH 2 density with a consequent increase in the 
axial velocity of the fluid. This effect then induces a large exit Mach number (M 
approaching 0.60 to 0.80) which shows that the bearing operates close to sonic 
conditions and may be accompanied by fluid vaporization. Dynamic operation 
of a liquid bearing at such large Mach numbers is entirely unknown, and thus, 
there is the urgent need to determinate experimentally whether test measurements 
produce a similar dynamic force response as the numerical predictions indicate. 

The barotropic-isothermal fluid model for turbulent flow HJBs has been 
extended to include journal misalignment and dynamic journal axis conical 
motions. The complex analysis determines rotordynamic force and moment 
coefficients due to journal center displacements and journal axis rotations (San 
Andres, 1993a). An approximate analysis for the HJB force and moment dynamic 
response at the centered position is also available (San Andres, 1992b). Numerical 
predictions show that journal misalignment angles in line with the journal center 
displacement cause a reduction in load capacity due to the loss in available 
film thickness. Journal axis misalignment increases slightly the flow rate and 
produces significant restoring couples for large angular displacements as well as 
large moment coefficients for large angular motions. The whirl frequency ratio 
for conical motions is equal to 0.50 in the absence of rotor gyroscopic effects, and 
therefore, identical to that of plain journal bearings. 

It is well known in fluid film lubrication theory and rotordynamics that there 
exists a trade off between the magnitude of the bearing static load capacity and 
its dynamic stability characteristics. In general, a bearing designer ’’sees” the 
need to reduce the former to improve the latter, or viceversa. Approximate 
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formulae for the calculation of a break frequency above which stiffness hardening 
and reduced damping occur in a compressible liquid HJB are available in the 
literature (San Andres, 1990b, 1991a, 1992b). Criteria to reduce the posibility 
of ’’pneumatic hammer” and to eliminate the reduced damping factor due to 
liquid compressibility effects are also formulated. In brief, improved dynamic 
characteristics (low values of the whirl frequency ratio) are attained for bearings 
with small recess volumes and very low values of the ratio between the pressure 
differential across the bearing and the liquid bulk modulus. Other recommended 
fixes to improve the dynamic characteristics of a hydrostatic bearing are: 

o Use of large scale roughened bearing surfaces to reduce the cross-coupled 
stiffness coefficients which are the direct leading agents in promoting hydrody- 
namic bearing instability. Experimental results at Texas A&M show that a test 
HJB with a rough knurled-pattern bearing surface has a WFR as low as 0.30 but 
with a reduced load capacity and direct stiffness when compared to a smooth 
surfaces HJB (Franchek et al., 1993). 

o Use of wear rings (Scharrer et al., 1992b) or end seal restrictions (San 
Andres, 1992b) to control bearing leakage, reduce the pressure differential across 
the bearing and increase the damping coefficients. Wear rings can also add a 
degree of safety for start-up and shut-down transient operation with potential shaft 
to bearing contact. The operating principle of end seals has been taken from the 
theory of shrouded pump-impellers (Childs, 1989). However, these wear ring seals 
in HJBs can lead to a substantial reduction of the direct force stiffness coefficient 
and static load support. Nonetheless, end seals of variable (controlable) restriction 
may prove in the future to be very valuable as active devices capable of producing 
load support and damping characteristics as desired within the environment of a 
’’smart” control rotor-bearing system. 
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o Use of tangential liquid injection opposing journal rotation to reduce the 
development of the circumferential flow velocity and to practically eliminate the 
cross-coupled stiffness coefficients. This concept yet lacks firm theoretical modeling 
but has proved to be useful in some applications while disastrous in others. 
However, recent experimental results at the hydrostatic bearing test facility at 
Texas A&M University have shown succesfully that a tangential injection HJB 
has a very low whirl frequency ratio without reducing the hydrostatic stiffness and 
load capacity (Franchek et al., 1993). These experiments add to our knowledge a 
conceptually simple fluid film bearing free of stability problems. Issues of concern 
for this especial bearing geometry are related to increased power losses and the 
posibility of inducing backward whirl in lightly loaded rotating structures. Braun 
et al. (1993) and Hill et al. (1993) have recently presented simplified numerical 
analysis relevant to non-radial fluid injection in laminar flow hydrostatic recess 
regions. 

o Bearing circumferential asymmetry has long been known to provide a 
minimal reduction in the whirl frequency ratio of cylindrical journal bearings. 
Elliptical and multi-lobe fixed pad bearing geometries have been the subject of 
past interest until the advent of tilt-pad bearings which are inherently stable. 
These bearings however are of large size, mechanically complex and present moving 
pivoted pads which need to be frictionless to operate efficiently. Furthermore, tilt- 
pad bearings present a dynamic response variable with frequency and coupled to 
the pads motion. So far these considerations have prevented conventional tilt- 
pad bearings to become suitable candidates as support elements in cryogenic 
turbomachinery. The flexure pivot tilt pad bearing introduced by KMC Inc. may 
overcome these problems (Zeidan, 1992). However, this type of bearing has yet to 
be tested for identification of its dynamic force coefficients. 
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Experimental results for the orifice discharge flow coefficients in a hydrostatic 
bearing (no journal rotation) for laminar and turbulent flows in freon, oil and 
gaseous nitrogen are given by Scharrer and Hibbs (1990). The experimental 
results point out to the complex nature of the flow field in the orifice and recess 
of a hydrostatic bearing and show that the orifice discharge coefficient is highly 
dependent on the Reynolds number with no evident correlation to the bearing 
clearance. The data presented is of utmost importance for the correct sizing of 
orifices for hydrostatic bearings in a cryogenic application. Scharrer, Hecht and 
Hibbs (1991a) present numerical results on the effects of bearing surface wear on 
the rotordynamic force coefficients in a HJB. The analysis is based on the solution 
of the turbulent flow Reynolds equation as given by the simplified model of Artiles 
et al. (1982). The results show that bearing performance is relatively unaffected 
by amounts of wear less than 10 percent of the bearing clearance. Direct stiffness 
and damping coefficients are found to decrease dramatically with increased wear 
while the cross-coupled stiffness is relatively unaffected. No results are presented 
in terms of the whirl frequency ratio to determine if bearing surface wear degrades 
the dynamic stability characteristics of a hydrostatic bearing. 

Scharrer and Henderson (1992a) present a detailed study on the specification 
and design of hydrostatic bearings for the STME fuel turbopump. The paper 
describes the important characteristics of the turbopump and addresses with 
detail to the major advantages of HJBs when compared to conventional roller 
bearings. HJBs are selected to provide maximum life and design flexibility at 
the minimum turbopump cost. A similar study for the design of a annular 
hydrostatic bearing to be incorporated into the SSME HPOTP is given also 
by Scharrer et al.( 1992c). This bearing package designed to act both as a load 
support element and a wear ring seal will replace the duplex ball bearing which 
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has been the primary life limiting component in the turbopump. The analyses 
of the static and rotordynamic characteristics of the hydrostatic bearings and the 
annular hydrostatic bearings are based on the numerical models developed by San 
Andres (1991c) and Yang et al.(1993a). 

Fundamental transient operation tests for liftoff and touch-down operation 
in a liquid nitrogen HJB are given by Scharrer, Tellier and Hibbs (1991b). A 
test apparatus simulating the transients found in a typical cryogenic turbopump 
was used in the investigation. Profiles of speed, pressure and radial load in the 
tester reproduce those determined from an instrumented SSME HPOTP. The 
experimental results show that a soft bearing material promotes early liftoff when 
used in conjuntion with a journal of hard coating. Post experimental inspection 
of the test especimens revealed unmeasurable wear for combinations of hard 
journal with soft bearing materials; and, conclusively that hydrostatic bearings are 
capable of surviving (repeated) start transients with similar loads as those found 
in actual turbopumps. A similar experimental procedure for the start transient 
performance of an annular hydrostatic bearing in liquid oxygen has been reported 
recently by Scharrer et al. (1992d). The experimental results showed that the 
fluid film bearing supported effectively the radial loads impossed and that liftoff 
occured at speeds very close to predictions obtained using the analysis of Yang et 
al. (1993a, b). 

Although cryogenic liquids offer very small viscosities, the trends toward 
higher rotational speeds and larger pressure differentials, as well as the implemen- 
tation of intentionally roughened surfaces to improve fluid film bearing dynamic 
stability, provide unique flow characteristics and operating conditions where high 
levels of turbulence (energy dissipation) and kinetic energy exchange due to fluid 
compressibility may determine significant thermal effects. Furthermore, hydro- 
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static annular pressure seals are also under consideration as support elements for 
cryogenic turbopumps. These type of damping bearings (Von Prageneu, 1982, 
1990) require particularly rough stator (bearing) surfaces and increased axial 
lengths to provide the necessary load support. These conditions then lead to the 
need of a thermohydrodynamic (THD) analysis in order to predict accurately the 
static and dynamic force performance characteristics of highly turbulent cryogenic 
liquid HJBs and seals. 

As stated earlier, cryogenic liquids present material properties strongly 
dependent on their local state of pressure and temperature. In some cryogenic 
turbopump applications, annular hydrostatic annular seals may be subjected to 
axial pressure drops as large as 45MPa, and therefore, the working fluid presents 
a wide spectrum of property changes across the seal. Pressure variations are 
most important on the viscosity and density changes of liquid hydrogen, while 
temperature has a pronounced influence on the material properties of liquid 
oxygen. The viscosity of liquid oxygen is an order of magnitude larger than 
in liquid hydrogen, but the kinematic viscosity is about the same for both fluids 
for the range of conditions applicable in cryogenic environments. Mechanical 
energy dissipation in turbulent flows is a non-linear increasing function of the fluid 
viscosity and the flow Reynolds number; and thus, thermal effects in liquid oxygen 
seals are expected to be of importance. This consideration becomes fundamental 
when it is further known that the heat capacity of liquid oxygen is relatively small. 

The understanding of thermohydrodynamic (THD) phenomena in hydrody- 
namic bearings has advanced steadily in recent years. A well documented review 
of the related literature relevant to THD lubrication theory and experiments is 
given by Khonsari (1987). The monumental work of Pinkus (1990) offers a com- 
prehensive review and insight of the current problems and describes the needs for 
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future work. The interested reader should refer to these publications for complete 
information. The conventional THD theory of hydrodynamic bearings consid- 
ers incompressible fluids of large viscosity (like oils), and thus, shows dominance 
of circumferential flow effects. Most analyses consider the flow as laminar and 
inertialess, and the axial variation of temperature is usually neglected. On the 
other hand, the THD analysis of externally pressurized bearings such as HJBs and 
damper seals needs to account for flow turbulence, fluid inertia and the transport 
of energy in the direction of the pressure flow. 

Isothermal- variable properties bulk-flow fluid film models for the analysis 
of annular seals in cryogenic environments have been adopted by Simon and 
Frene (1989) and San Andres (1991c, 1993b-c). These models still account for 
the material properties variation and have been shown to be accurate for design 
of liquid hydrogen seals as described latter. Yang et al. (1993a, b) have presented 
the most advanced THD analysis for the turbulent bulk-flow of a real properties 
liquid in an annular pressure seal of arbitrary non-uniform clearance. The thermal 
model considers conservation of bulk-flow momentum and energy transport while 
an adiabatic heat flow process is assumed for simplicity. The bulk-flow equations 
for the zeroth-order flow variables are solved numerically to determine leakage, 
fluid film forces and restoring moments. Dynamic force and moment coefficients 
due to small amplitude perturbations in shaft center displacements and shaft 
axis rotations are calculated from the solution of first-order flow field equations. 
The material properties used for cryogenic liquids are determined from the 32- 
term Benedict-Webb-Rubin equation of state as given in the standard computer 
program of McCarty (1986). 

Numerical results from the THD model correlate well with measured temper- 
atures for the Preburner Impeller Rear Wear-Ring Seal of the SSME high pres- 
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sure oxidizer turbopump. Calculations performed with the analysis of Yang et al. 
(1993a) for a high speed, rough surface LO 2 seal predicted a large temperature 
rise and the onset of two-phase flow conditions (liquid boiling) at moderate shaft 
eccentricities despite the large pressure differentials impossed across the seal. It 
is imperative to note that large temperature rises in a cryogenic seal can lead to 
thermal solid distortions affecting the operating clearance and possibly causing a 
significant reduction on the direct stiffness of the hydrostatic seal. 

Two-phase (liquid-gas) flows are not desirable in cryogenic liquid fluid film 
bearings and seals. A quote from Pinkus (1990) is the more explicit since ’’difficul- 
ties attendant to the phenomenon of change of phase includes a negative stiffness 
over certain ranges of operation, leading to unstable equilibrium and collapse of 
the fluid film. Even when operation is possible, self-sustained oscillations may lead 
to chatter, ultimately causing destruction of the seal”. Pinkus also states that the 
analysis of two-phase lubrication with cryogenic liquids is extremely complex due 
to the ’’inadequate state of knowledge of certain branches of fluid dynamics in 
general and the lack of basic experiments”. Beatty and Hughes (1987, 1990) have 
introduced fundamental studies on two-phase flow phenomena for steady state, 
centered operation, annular seals with adiabatic surface considerations. The anal- 
yses showed that viscous heat generation is the mechanism for limiting seal leakage 
as the operating conditions are near two-phase regions. 

Research is needed to develop a thermohydrodynamic model for the analysis 
of turbulent single-phase and two-phase flows in hydrostatic bearings and annular 
seals for cryogenic applications. Advanced analytical and numerical tools are 
needed to predict accurately the steay-state and dynamic force response of rough 
surface, fluid film bearings with significant levels of energy dissipation. Typical 
hydrostatic bearings and annular seals in a cryogenic environment operate with 
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large pressure differentials and, in some circumstances, the inlet temperature to 
the bearing/seal is just a few degrees below the liquid saturation temperature. For 
example, hydrostatic annular seals control effectively leakage and could also (in 
the near future) replace conventional rolling element bearings as load support 
devices in high speed turbopumps. Thus, the design of these seals (damper 
bearings) calls for operation at moderate to large eccentric rotor positions in 
order to generate adequate load support. On the other hand, these operating 
characteristics generate the conversion to heat of mechanical energy due to the 
hydrodynamic effects as well as that due to the macroscopic roughness of the stator 
surfaces used. Recent experimental evidence from tests beds at both Rockwell and 
NASA Marshall Research Center (Nolan et al., 1993) have shown that prototype 
damper seals present a two-phase flow condition (flashing) of the liquid at the 
discharge plane of the damper seal. Operation of actual turbopumps under these 
conditions may be potentially catastrophic. 
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SCOPE OF THE WORK IN 1993 (PHASE I) 

The research progress presented in this report concentrates on the thermohy- 
drodynamic analysis of real properties, fluid film bearings for cryogenic applica- 
tions. The bearing geometries of interest include hydrostatic bearings of arbitrary 
recess geometry, hydrostatic annular seals and cylindrical rigid-pad bearings. The 
motion of a cryogenic fluid on the thin film flow region is governed by a closed 
set of mass and momentum conservation, and energy transport equations for the 
primitive bulk-flow variables and, accompanied by thermophysical state equations 
for evaluation of the fluid material properties. The turbulent bulk-flow in the thin 
film land regions is modeled with simple friction coefficient laws including the 
effect of macroscopic surface roughness. Zeroth-order equations describe the fluid 
flow field for a journal static equilibrium position, while first-order linear equa- 
tions govern the fluid flow for small amplitude journal center translational motions 
and journal axis conical motions. Solution to the zeroth-order flow field equations 
provides the bearing flow rate, load capacity, restoring moments and torque, and 
the temperature field within the flow field. Solution to the first-order equations 
determines the linearized force and moment coefficients due to journal lateral and 
angular motions. 

On 1993 (Phase I), the analysis of single-row and side-to-side parallel row, 
rectangular recess hydrostatic bearings, annular pressure seals, and fixed geometry 
journal pad bearings has been completed. Thermal effects (energy transport) 
are considered in the fluid film bearing with thermally adiabatic solid surfaces 
bounding the annular flow region, or isothermal journal and bearing surfaces. 

The computer program developed named as hydrosealt calculates: 

1) bearing flowrate (seed leakage), 

2) friction torque and temperature rise, 

3) load capacity (fluid film forces) and restoring moments, 



18 


4) rotordynamic force coefficients due to journal center displacements, 

5) rotordynamic force coefficients due to journal axis rotations, 

6) rotordynamic moment coefficients due to journal center displacements, 

7) rotordynamic moment coefficients due to journal axis rotations, 

8) Complete pressure and temperature fields on the bearing surface, as well 

as density and viscosity field variations, with ranges of fluid flow 

Reynolds numbers and Mach numbers. 

for isothermal or adiabatic journal and bearing surface conditions. 

as a function of: 

a) journal (rotor) center eccentricity and journal axis misalignment. 

b) inlet specified circumferential pre-swirl velocity distribution. 

c) general clearance function as defined below 

d) mean surface roughness on bearing and journal film lands. 

The fluid properties (density, viscosity and specific heat) are calculated from 
the 32 term Benedict- Webb-Rubin equations of state as given by NBS Standard 
Reference Database 12 (McCarty, 1986) for the following cryogenic liquids: LO 2 , 
LH 2 , LN 2 , and methane. Other fluids included in the program are water, mineral 
oils and air. 

The program hydrosealt handles the following boundary conditions at the 
bearing discharge (exit) planes: 

1) periodic pressure asymmetry in the circumferential direction. 

2) local discharge end seal effects via an orifice like model to simulate 

end wear-rings in hydrostatic bearings or annular seals. 

The axial clearance functions included are of the type: 

a) uniform, 

b) tapered, 

c) stepped, or, 

d) arbitrary via spline interpolation, and, 

e) tilted about any arbitrary axis for journal misalignment conditions. 

The present computational research program needs support from a compli- 
mentary experimental research effort directed to measure the dynamic force per- 
formance of high speed, externally pressurized fluid film bearings with flow con- 
ditions similar to those found in cryogenic environments. Some empirical flow 
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parameters should also be identified from the test data. These coefficients, like 
orifice discharge coefficients, recess edge entrance loss factors, shear friction coeffi- 
cients, and thermal mixing groove parameters are of fundamental importance for 
the accurate prediction of the flow and force dynamics in high-speed, turbulent 
flow fluid film bearings. 



20 


ANALYSIS 

Fluid Film Bearing Geometries 

The following bearing geometries are considered in the analysis: 

(1) Orifice-compensated hydrostatic journal /pad bearings (Fig.l). 

The type of bearing being considered for use in the LHi and LOi turbopumps 
is a full, 360° hydrostatic journal bearing, orifice compensated, with a variable 
number of feeding recesses or pockets machined in the surface of the bearing. 
Axial grooves may be designed in the bearing surface to improve the bearing 
dynamic stability. The bearing surface may also be intentionally roughened to 
reduce the cross-coupled stiffnesses and improve the bearing stability. 

The bearing recesses have a depth larger than the radial clearance under the 
film lands. Fluid flows from a high-pressure supply, passes through an orifice in 
route to the recess and leaves the bearing sides by flowing through the thin film 
gap. The surfaces of the journal and bushing are separated by a film of fluid 
forced under external pressurization. When the center of the shaft is displaced 
from the center of the bearing, on one side of the bearing the film gap is reduced, 
while the gap is increased on the other side. On the reduced-gap side, the flow 
from the recess to film lands is more restricted and the recess pressure rises. On 
the opposite side, however, the flow from the recess is increased and the recess 
pressure decreases. The difference between the pressures on the two sides results 
in a net bearing film force which is equal and opposite to the externally applied 
load on the journal. The hydrodynamic effect due to journal rotation increases the 
fluid film forces but also introduces the potential of a self-excited whirl instability 
at rotor speeds twice above the system first critical speed. 

(2) Annular pressure (damperl seals or hydrostatic annular bearings (Fig. 2). 

Annular pressure seals are used in turbomachinery to limit the leakage 




Fig.l Geometry of Hydrostatic Pad Bearing 
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of fluid from stages of high pressure to those of lower pressure. Like fluid 
film journal bearings, annular pressure seals have a significant influence on 
the rotordynamic stability of high-performance turbomachinery. Von Pragenau 
(1982, 1990) proposed damping seals, i.e., annular seals using rough stators and 
smooth rotors, to improve rotordynamic stability. In a damping seal, the average 
tangential velocity induced by journal rotation and the cross-coupled stiffnesses 
decrease. Scharrer et al.( 1992a) showed that a damper seal (called a hydrostatic 
annular bearing there) can also be used as a supporting element to replace a 
current rolling-element bearing in cryogenic turbopumps. 

(3) Cylindrical pad journal bearings (Fig. 3). 

The working mechanism of a hydrodynamic journal bearing is different from 
a HJB. The pressure generation in a hydrodynamic bearing is provided by journal 
rotation and fluid viscosity. External pressurization serves only to provide enough 
lubricant to the fluid film region. Therefore, a hydrodynamic journal bearing is 
also called a self-acting bearing. Unlike a hydrostatic bearing, a hydrodynamic 
bearing has no load capacity unless the journal rotates at an off-center position 
(a wedge like hydrodynamic effect). The lack of hydrodynamics may lead to wear 
damage during startup and shutdown transients. 

For fixed pad journal bearings, feeding axial grooves are located at circum- 
ferential locations away from the maximum pressure (minimum film thickness) 
region. Feeding grooves supply fresh (cold) fluid to the bearing, i.e., provide the 
required lubricant flow to the fluid film region and are critical in avoiding journal- 
bearing overheating at high speed operation. Due to the low level of external 
pressurization (relative to ambient conditions), the mass flow rate in a hydrody- 
namic journal bearing is much smaller than in a HJB or an annular pressure seal, 
and consequently thermal effects may be of great importance for fluids with low 
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heat capacity. A too large temperature rise in the fluid film affects pronouncedly 
the fluid properties, and hence the bearing performance characteristics. More im- 
portantly, unequal thermal growth of bearing components can modify the fluid 
film clearance and induce seizure. Bearings with compliant surfaces (foil bear- 
ings) exhibit advantages over rigid surface bearings by accommodating thermal 
distortion or mechanical deformation of bearing components and journal misalign- 
ments. NASA currently sponsors the development of fluid-film foil bearings for 
use in space shuttles and military missiles (O’Connor, 1993 ). 

For all the three fluid film bearing configurations considered above, the flow 
is confined to the thin film annular region between an inner rotating journal and a 
non-rotating bushing. The fluid flow is characterized by high levels of turbulence 
due to the high journal surface speed and/or the externally imposed large axial 
pressure drop across the bearing/seal. Two different flow patterns are defined in 
the analysis. 

Flow On Film Lands 

Fully developed, single-phase (liquid or gas), turbulent bulk-flow model; 

Governing Equations: 

Mass conservation equation 
Momentum Bulk-Flow equations 
Energy Transport equation 

Turbulence Closure Model: 

Bulk flow with friction parameters based on Moody’s 
friction-factor equation for physically sound characterization 
of roughened surfaces. 

Liquid of variable properties as functions of pressure and temperature with 
realistic thermophysical equations of state 
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(for cryogenic fluids, see McCarty, 1986). 

Flow On Bearing Recesses and Grooves 

Global mass conservation equation with liquid compressibility at recess- volumes. 
Pressure rise/ drop at edge of bearing recesses with a simplified 
1-D Rayleigh step bearing model. 

Global energy balance equations at bearing recesses and grooves to determine 
entrance temperature conditions at the recess and boundaries. 

Governing Equations for Turbulent Bulk-Flows 


The motion of a fluid in thin film geometries is described by a set of mass, 
momentum conservation, and energy transport equations for the primitive flow 
variables, and accompanied by thermophysical state equations for evaluation of 
the fluid material properties. The smallness of the clearance to radius ratio (c*/i2) 
allows the reduction of the general governing equations to the simplified forms 
given on the film lands (\ang et al., 1993a): (Note: the nomenclature of this 
report is provided in Appendix A) 

Continuity Equation 


dp dpU dpV dpW 

Tt + IT + + ~dT = 0 

Circumferential- Momentum. Equation 

DU = _dP dr zz 

^ Dt dx dz 


( 1 ) 


( 2 ) 


Axial- Momentum Equation 


DV 

'~Dt 


dP dr , 


dy 


+ 


yz 


dz 


( 3 ) 
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Cross- Film- Momentum Equation 


dP 

dz 


= 0 


Energy- Transport Equation (Bird et al.,1960 , Whitaker, 1977) 




( 4 ) 


( 5 ) 


on the flow region {x E (0,rrD),y E ( — Ll,Lr),z E (0, H(x, y, <)), for journal 
bearings and hydrostatic bearings, while y E (0, L) for annular pressure seals}; 
and where 


D 

Dt 




( 6 ) 


is the material derivative. 

To undertake a study of the general equations (1) — (5) in their full three- 
dimensional complexity is a difficult task, perhaps unnecessary to obtain mean- 
ingful solutions in a reasonable time and at an affordable cost. Large pressure 
gradients, typical of externally pressurized cryogenic bearings and seals, and high 
rotational speeds generate large flow Reynolds numbers, and consequently, the 
effect of turbulent mixing far outweighs molecular diffusivity. Thus, the tempera- 
ture rise produced by viscous dissipation tends to be distributed uniformly across 
the film thickness and hence temperature gradients in the cross-film coordinate 
(z) are confined to turbulent flow boundary layers adjacent to the bounding (bear- 
ing and journal) surfaces (Suganami and Szeri, 1979; Di Pasquantonio and Sala, 
1984). Furthermore, in the absence of regions of reversed flow or recirculation, 
the fluid velocity field presents the same characteristics as discussed above. 

Let bulk-flow primitive variables be defined as average quantities across the 



film thickness, i.e. 
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1 f H - 1 f H . 1 f H . 

u= «L uiz ’ v= nL vdz ' t= hJ 0 Tit < 7 > 

Then, the considerations presented allow the governing equations (1) to (5) to 
be integrated across the film thickness to determine turbulent bulk-flow equations 
at each bearing pad configuration (Yang et al., 1993a): 

Continuity Equation 


d(pH) djpHU) d(pHV) 

dt dx dy 


Circumferential- Momentum Equation 


djpHU) _ djpHlP) djpHUV) 
dt + dx + dy 


rjdP 

~H T Txz 

ox 


H 

0 


Axial- Momentum Equation 


(8) 


(9) 


d(pHV) d(pHUV) d{ P HV 2 ) TT dP 

dt + dx + dy “ ~ H ~dy + T ' 


H 

yz lo 


( 10 ) 


Energy- Transport Equation 


r d(pHT) d(pHUT) d{pHVT)i 


+ 


+ 


+ 2L 


dt dx ' dy J 

lr d(pHV?) d(pHUV?) , d(pHVV ?) i 


dt 

+ (T/ 3, - 1 )H(U^ + V^)+ 


+ 


dx 


+ 


dy 


+ Qt 


(ii) 


Note that Eqs.(9) to (11) imply “flat” bulk flow profiles characteristic of 
turbulent film flows at high Reynolds numbers. By Neglecting the fluid-inertia 
terms on the left hand sides of Eqs.(9) and (10), Eqs.(8) to (10) reduce to a single 
lubrication equation — namely, the turbulent flow Reynolds equation. 
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In the above equations, V t = y/U* + V 2 is the bulk-flow speed, r xz |g 
and r yz iff are the wall shear stress differences in the circumferential and axial 
directions, t xz \ h is the shear stress on the journal surface in the circumferential 
direction, and 

Q* = h g (T -T b ) + hj(T - T s ) (12) 

is the heat flux from the film to the bounding (bearing and journal) surfaces. 
Further, h Q and h 3 denote the bulk-flow heat transfer coefficients to the bearing 
and journal surfaces (Appendix B). If the fluid-film temperature ( T ) is higher 
than the bounding-surface temperatures (T b ,Tj), heat flows from the film to the 
bounding solids, and vice versa. For an adiabatic flow process, Q s = 0. 

The fluid properties for a cryogenic liquid (density, viscosity, specific heat, 
and volumetric expansion coefficient, etc.): 

p = p(P,T), p — p(P,T), C p = C p (P,T), (3 t =(3t(P,T) (13) 

are determined from general 32-term thermophysical state equations as given by 
McCarty (1986). 

The wall shear stresses are calculated according to the bulk-flow theory for 
turbulence in thin film flows (Hirs, 1973, Launder and Leschziner, 1978): 

r yzlo (14) 

Txz ^ H = + W [Uka ~^ U ~ 

where the turbulent shear parameters (k x ,k y ) and (k ,k B ) are local functions of 
the Reynolds numbers and friction factors based on Moody’s formulae (Massey, 
1983). 



30 


Substitution of the bulk-flow momentum equations (9), (10) into the energy- 
equation (11), and invoking the mass conservation principle given by equation (8) 
allows the energy transport equation to be expressed as: 

n ApHT) d(pHUT) , d(pffVT) , „ 




n T |» _ Vt I" 

v T xz lo v T yz io 


which reflects the energy balance in the fluid film as: 

CONVECTION + DIFFUSION = COMPRESSION WORK + DISSIPATION 


( Energy Disposed ) 


( Energy Generated 


The variation of temperature in the axial direction and the energy generated 
by compression work are retained in the analysis due to the strong influence of the 
large pressure drop across the length of a HJB or an annular pressure seal. These 
conditions differentiate the present problem from conventional THD analyses of 
viscous, incompressible fluids in journal bearings where the volumetric expansion 
coefficient (3 t is zero and temperature variations along the bearing axial direction 
are usually considered negligible. 


Governing Equations for Bearing Recess and Groove Flows 


The analysis of turbulent flows in a HJB recess is complicated and not yet 
fully understood. To date, only two-dimensional laminar flow numerical solutions 
are available for rectangular recesses (See, for example, San Andres and Velthuis, 
1992, Hill et al., 1993, Braun et al. 1993). While the actual prediction of flow 
fields in the recess may give a better description of the recess-film edge boundary 
conditions, the global mass and energy conservation principles at the recess are 
shown to be both efficient and accurate in HJB analysis. The study of the flow field 
in a bearing recess is a subject which requires detailed experimental results since 
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analytical or full 3-D numerical solutions can not model efficiently the complexity 
of this flow field dominated by both surface shear (rotation) and induced external 
pressure flows of strong turbulence character. 

• Mass Conservation at a Recess 

The continuity equation at the recess is defined by the globed balance between 
the supply mass flow rate (Q; n ), the outflow into the film lands ((? r ) and the 
temporal change of fluid mass within the recess volume (V r ). The recess flow 
continuity equation is expressed as: 

Qin = Qr + P r^+ p,V r (/? p ^ - Pt~) r ( 16 ) 

where the flow through the orifice restrictor is 

Qin = A 0 vV(P 3 - P r ), (17) 


and 


Qr= f pH(\J ■ n) dr (18) 

J r, 

is the mass flow rate across the recess boundary (r r ) and entering the film lands. 



0t 


p\dT / p 


(19) 


are the liquid compressibility factor and volumetric expansion coefficient, respec- 
tively. In general, 


{ = 0 for incompressible liquids; 

= p, 5 ; for ideal gases; ( 20 ) 

> 0 for cryogenic liquids. 

In HJBs, the compressibility of cryogenic liquids has been shown to promote 
pneumatic hammer, a class of self-excited instability mechanism as described by 
San Andres (1991). 
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• Global Energy Balance Equation at a Recess 

The energy transport phenomenon in a HJB recess is controlled by the follow- 
ing three mechanisms: the carry-over of hot fluid from upstream to downstream of 
the recess, the mixing of fresh fluid from the supply source into the recess volume, 
and the heat generation in the recess volume due to shear dissipation by journal ro- 
tation. Energy transport produced by pressure gradients, kinetic energy changes, 
and heat conduction are considered negligible due to the uniformity of pressure 
and the large mass flow rate through the recess. Based on these considerations, 
the general three-dimensional energy equation (5) reduces to 



DT d - 

pC - D t - 9^ Ut » + Vt »■> 

(21) 

Invoking the continuity equation (1), Eq.(21) can be written as 


c P 

[£ +v-(„vr)] _ £-(Ut„ + Vt, z ) 

(22) 

where 

V = {U,V,W} 

(23) 


Integration of Eq.(22) over the recess volume gives the global energy balance 
equation which reflects the heat carry-over (advection) and mixing effects, and 
the friction heat generation (dissipation) in the recess (Fig. 4) (C p is regarded as 
constant for simplicity in this part of the analysis): 


Cp ^F Vr + Cp (S rhdTd + £ rh 3 ide T 3ide ) 
= c p (J2 ™uT u + Q in r ,) + r; r n 


where 


(24) 


™'d,u — (pU HSy)^^ and = (pV H8x) sic i e (25) 

are the mass flow rates over edge differential segments. The summation (XI) is 
over each recess edge, and the subscripts “u”, “d” and “side” refer to the upstream, 
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downstream, and side (axial) edges of a rectangular recess, respectively, 
the total mass flow rate through the supply orifice, and 

T„ = r»A r R 


Q in Is 


(26) 


is the drag torque over the recess with volume V r . The left hand side of Eq.(24) 
represents the temporal variation of energy in the recess volume and the energy 
flow rate out of the recess, while the right hand side denotes the energy influx and 
the dissipation of mechanical energy. 

Our purpose here is to obtain the recess-edge temperatures which serve as 
boundary conditions for solution of the thermal field on the film lands. The 
temperatures at all the recess edges and within the recess can be obtained from 
the global energy equation for different recess flow conditions. Since the global 
energy equation involves the time (i) as a parameter, different treatments are 
needed for the zeroth- and first-order temperatures produced by the perturbation 
analysis. Appendix C provides the formulas for calculation of the zeroth- and 
first-order recess/edge temperatures. 

• Global Energy/Flow Balance Equation at a Groove 

Axial bearing grooves (Figs.l or 3) are used for different purposes, such as to 
feed the bearing, to improve stability of a HJB by providing bearing asymmetry, 
or to remove foreign particles in the lubricant. A groove is usually narrow (small 
circumferential width) and deep (depth much larger than the film thickness). 
Therefore, the grooves are treated separately from the film lands. 

At an axial groove, the flow mixes and energy transfer occurs with cooling 
of the pad leading edge flow if the bearing is properly designed. The problem is 
complicated and subject to empiricism. Here a global model using averaged flow 
variables along the groove edges is implemented. 



Journal Surface 



(b) Qi n >0 (C) Q i n <°' Q a >0 (d) Q in <0 ' Q d <0 


Fig. 5 Conceptual description of global energy 
balance at a groove 
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Fig. 5 presents the conceptual description of global energy balance at a 
groove. The energy transport phenomenon in a groove is controlled by the carry- 
over of hot fluid from the upstream pad to the downstream one, and the mixing of 
fresh fluid from the supply source into the groove. The global energy conservation 
equation at a groove is given as ( C p is regarded as constant here): 

C P ^™ d T d = C p ^rm u T u + C p Q in T in (27) 

where m d and m u are the mass flow rates over edge differential segments as 
given by Eq.(25), and the summation (V)) is over the downstream (leading) and 
upstream (trailing) groove edges. Also, 

Qin = dip? - y~^ m u — Q d ~ Qu (28) 

represents the net mass flow rate through the groove, providing the minimum 
flow rate required to prevent starved lubrication for a journal bearing with a 
single feeding groove. T d and T u represent the temperatures at the downstream 
and upstream groove edges, and T; n is the temperature flowing with mass flow 
Qin- 

If the mass flow rate leaving the trailing edge of the upstream (relative to the 
journal rotation) pad to the groove is greater than zero, i.e., 

Qu - ^ rh u > 0, (29) 

then, the fluid temperature at the trailing edge can be found by their upstream 
values given by the solution on film lands (upstream pad): 

T u = Upstream Film Values on Upstream Pad (30) 

Temperatures at the downstream edge (leading edge of the downstream pad) can 
be obtained for different flow conditions in the groove: 
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(a) If Qj = E ril'd ^ Qm Q in ^ 0 and flow is supplied into the groove, 
then the flow entering the downstream edge of the groove is the combination of 
the upsteam (Q u ) and the groove supply (Qi n ) flows (Fig. 5b), and an averaged 
uniform downstream-edge temperature (T^) can be obtained from the mixing of 
upstream and inlet temperatures: 

\Vrii u T u / Qu\ 

= + (i-a|)t,„ ,3i) 

where 0.4 < A < 1.0 is the mixing coefficient, a participation ratio of the hot 
lubricant (recirculating fluid) in the mixing process (Mitsui et al., 1983). 

Tin = T a (32) 

that is, all make-up flow in the groove is at the fluid supply temperature. 

(b) If Qd < Qu? i.e., Q ‘s 0 and fl ow is discharged out of the groove, then 
two cases may occur. First, a fraction of the upstream flow enters the downstream 
pad ( Qd > 0)(Fig.5c). Then, the leading edge temperatures(Td) are equal to the 
trailing edge ones ( T u ). Second, both edge flows enter the groove ( Q u > 0 and 
Qd < 0) (Fig.5d) and the temperatures on both edges can be obtained from 
thermal solutions on their respective pads. These two cases are typical of HJBs 
and represented by: 

T d ={ ' KTu ’ if Qd > o (Fig. 5c); . 

( Upstream Film Values on Downstream Pad, if Q d < 0 (Fig.5d). 

• Recess /Film Entrance Pressure Rise/Drop at HJBs 

For purely hydrostatic operations, a uniform recess pressure is desirable 
not only to simplify the HJB design (and analysis) but also to increase the 
load capacity. A uniform recess pressure can be achieved by deepening the 
recess. However, a minimum recess volume is required to avoid the characteristic 
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pneumatic hammer instability problem in HJBs. Design criteria for uniformity of 
recess pressure and pneumatic hammer instability are given by Redecliff and Vohr 
(1969), and San Andres (1991a, 1992b). 

For hybrid operation of a HJB, a pressure rise due to the journal rotation is 
produced in the downstream portion of the recess (Ho and Chen, 1984, Chaomleffel 
and Nicholas, 1986). San Andres (1992a) considered this region as a one- 
dimensional step bearing and evaluated the pressure rise just in front of the 
downstream recess edge (P~) as (Fig. 6 ): 


P 7 = 


where 


{ p - 

brf 2 y,k r 


{ 

2H^(l-iV/ r 2 ) 

{ Pr, 




in y direction 
(34) 


U — Ui + Vj, V r — RQi -|- 0 j, 


(35) 


- ^ 0 001 ./ * 

= R& r /7 .753 is a turbulent shear factor at the bearing recesses, 77 = 
H/(H r + H) is the ratio of thicknesses at film land and recess depth, and 
M r = OiRQy/ ]Tp~ r is the recess flow Mach number for hybrid operation. 

The local acceleration of fluid from the deep recess to the film lands causes a 
sudden pressure drop at the recess edges (Fig. 6 ). The pressure at the entrance to 
the film lands is modeled by simple Bernoulli type relations based on the turbulent 
flow theory developed by Constantinescu and Galetuse (1975): 


- §■(! +<«)(! + £»)[! - 7 2 ( — )](U 


n) 2 ; (U • n) > 0 


(36) 


(37) 


where 


(i = {CitCy} = {l-95/i2e°' 43 , 0}, 


( 38 ) 
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Fig. 6 Conceptual description of pressure rise 
and drop at recess edge of a hydrostatic 
pad bearing, and pressure ram effect at 
leading edge of bearing pad. 
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are the empirical entrance loss coefficients and n is the normal vector to the recess 
boundary (T r ). 

The pressure at the entrance of an annular pressure seal can be obtained 
similarly by using a Bernoulli type relation. The corresponding expression is 
given latter in the boundary-conditions section. 

Dimensionless Governing Equations at Film Lands 

Dimensionless coordinates and variables are defined as follows 


R 1 h= T.' T= “' ; 


_ V V p-p, 

U TT ’ V TT > P “ An 


p = ~~ ; p = — ; c P = p 

P* A 4 * C p 


- T 
T = — 

T* 


= PtT * ; 


and TT — c « I ±t_ 

U * ~ ti.R 

is a characteristic speed due to pressure induced flow with A P = P, - P a for 
externally pressurized bearings (HJBs and annular pressure seals). If P 3 — P a , 
like in a pure hydrodynamic bearing, the characteristic speed is set equal to the 
journal surface speed 

U , . = Rn MOV 


= vm(— 

\c- 


In dimensionless form, the governing flow equations at the film lands 


written as: 



( 42 ) 
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Circumferential- Momentum Equation 


dp 


d 


h dx~h^ kzU kj 2 ^ + Rea ^(P hu ) + Re pi^(P hu2 )+ -^(phuv)}-, (43) 


dr 


dx 


dy 


Axial- Momentum Equation 


h dy~ + Re ^(P hv ) + R^p[^r(phuv) + £=(phv 2 )}‘, 


Energy- Transport Equation 


C p Re » q t {phT) + C p Re*[—(phuT ) -f phvT )] + Re*Q 3 


l dx' r -—' ' dy 

— r? $ 4 u( dp dp x - A dp 

—E c \ dth(<T — — |- u— + v--)T + h — — — 
1 dr dx dy' 2 dx 

+ ~~[k x {v\ 4- - uA) + ^(-A 2 - uA)] j 
Xhe dimensionless flow parameters are defined as 


RQ U\ 

u. ’ u. ; c “ r.c p . 5 


D. _ D • D C, , 

"*•- = aRe r' Re >=-^r ! Re P = Re p(j) 


p* 


(44) 


(45) 


In externally pressurized bearings the frequency (<r) and speed (A) numbers 
denote the importance of squeeze film and hydrodynamic flow effects relative to the 
pressure induced flow, respectively. The reference Reynold numbers (Re p ) denotes 
the ratio of fluid advection forces to viscous flow induced forces. In most HJBs, 
flow turbulence is produced by the Poiseuille flow due to the high pressure drop 
across the bearing length. The ratio of ( Re p /E c ) or ( Re,/E c ) may be interpreted 
as the effect of heat convection relative to shear dissipation. The Eckert number 
(£ c ) denotes the ratio of kinetic energy to heat convection in the fluid film. 


Dimensionless Governing Equations at Recesses 
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In dimensionless form, the governing equations at recesses are written as: 


• Mass Conservation at a Recess 



- dv r _ 

<lin — <?r + &Pr Pr ( 

( 5 dp r 5 0T r \ 

A-ar A 'c>J 

(46) 

where 



9tn = ^+\/pr{ 1 Pr) 

for orifice flow 

( 47 ) 

q r = J ph( u • 

n) dT r , 

(48) 

6 , = A 0 fi » v2 ![c\ 

VpTap), 

(49) 

II 

<1 

‘ } 

(50) 

fipr — (3 pr AP , 

(51) 

fltr — fitrT* , 

(52) 

* 

O 

> 

II 

i> 

R 2 ) 

(53) 


The orifice parameter 6 m is a design variable to determine a certain recess 
pressure ratio at the desired operating condition, and thus determines the size of 
the orifice restrictor. The major uncertain factor in the orifice parameter 6» is 
the value of the orifice discharge coefficient Cd which is highly dependent on the 
local Reynolds number in the orifice feeding hole as well as on downstream flow 
conditions. In general, Cd takes values ranging from 0.6 for sharp edged orifice 
configurations and highly compressible fluids, to 1.0 for highly turbulent flows and 
reduced orifice length-to-diameter ratios (Redecliff and Vohr, 1969). However, 
the experimental measurements performed by Scharrer and Hibbs (1990) show 
that the discharge coefficients can be as low as 0.3 for large pressure drops in 
hydrostatic bearing geometries handling cryogenic fluids. Extensive experiments 
show that the discharge coefficient (Cd) takes values between 0.75 and 0.95 for 
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water HJBs with different clearance, orifice diameters, and operating conditions 
(Kurtin et al., 1993, Franchek and Childs., 1993, and Mosher, 1993). 

The frequency parameter at a recess volume (A a ), like cr, denotes the ratio 
of squeeze film flow to the pressure induced flow. /3 pr is the dimensionless recess 
volume-liquid compressibility factor, which should be small to avoid pneumatic 
hammer instability at low frequencies (Reddecliff and Vohr, 1969, San Andres, 
1991). A combination of large compressibility and frequency parameters, f3 pr and 
Aj, could bring undesirable hardening stiffness and substantial loss of damping 
at high excitation frequencies (Rohde and Ezzat, 1976, Ghosh et al., 1979,1987). 
Note that the introduction of temperature variation in the recess may mitigate or 
deteriorate the pneumatic instability problem depending on the sign of (see 
Eq.(46)). 


Global Energy Balance Equation at a Recess 

d(pT) 


rhdTd + ^ tilsideTside^ 

= + 9'") + (f; r A )EJR, 


(54) 


where 


rpr 

* or 

° r ~ c.APIP ' 


(55) 


is the dimensionless torque over the recess area. 


Global Energy Balance Equation at a Groove 


E = E ™-uT u + qinTin 


(56) 


where 


9in “ ~ Win — Qd — <lu 


(57) 



If flow leaves the traling edge of a pad then 
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<lu 


= £ 


m u > 0 , 


(58) 


so that 


T u — Upstream Film Values on Upstream Pad (59) 

and the dimensionless temperatures at the downstream edge of a feeding groove 
are given by: 

(a) If <ld = X) m d > q Uy then 


where 


T d 


b 

Qd 



Tin = 1 . 0 ; 


(•>) If«J < q u , then 


(60) 


(61) 


T d 


XT U or Ti n , 

Upstream Film Values on Downstream Pad, 


if q d > 0; 
if q d < 0. 


(62) 


Recess/Fi lm Entrance Pressure Rise/Drop 


p ~ - J Pr ~ Pr<t>r 
{ Pr, 


^(fy) ~ yj • n, (u • n) > 0, in x direction 

in y direction; 


(63) 


+ — 4>iPe / _> o / -> 

Pe =Pe r-(u.n) 2 ; (u • n) > 0 


(64) 


where 


u = ui + vj, A = A i + 0j, (65) 

4>r = k r (b/D)( v /h) 2 /(l-M 2 r ) (66) 

4>i = <f> Xty = Re' p ( 1 + C, t ,)(l + &,„)[! - ( VPe/pr ) 2 ] (67) 
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Perturbation Analysis 

A perturbation analysis separates the governing equations into zeroth- and 
first-order equations. The zeroth-order equations describe the fluid flow field 
for a journal static equilibrium position, while the first-order linear equations 
govern the fluid flow for small amplitude journal center translational motions and 
journal axis conical motions. Solution to the zeroth-order flow field equations 
provides the bearing flow rate, film forces and load capacity, restoring moments 
and torque. Solution to the first-order equations determines the rotor dynamic 
force and moment coefficients due to journal lateral and angular motions. 

Fig. 7 illustrates the journal translational and rotational motions in an 
annular seal. The inertial coordinate system ,Z} helps to define the position 

of the spinning journal. For steady state operating conditions, the journal is at 
an equilibrium position given by two journal center displacements (e ro , e^ 0 ), and 
by two angular rotations (6xo> <$yo) about {X,T} axes and centered at the axial 
coordinate Zo, as shown in Fig. 7. Note that the origin of the inertial coordinate 
system in a HJB is different from that of a seal. For HJBs, the Z axis starts at the 
mid- plane of bearing recesses (San Andres, 1993a). The angles S x o or 8y 0 define a 
misalignment condition for the journal axes relative to the stator. Superimposed 
on the static equilibrium position, the journal describes motions of small amplitude 
(Ae x , Ae^, Fig. 7a) as well as angular rotations of amplitude (A# x , A <5 y , Fig. 7b) 
with center at Z q. The film thickness is given by the superposition of steady-state 
(fio) and dynamic (/ii) components, and represented in dimensionless form by the 
real part of the following expression: 

h = ho + hi (68) 


ho — cq + [e XQ -f Syo(Z — Z 0 )/c,] cos# + [e Yo — 6xo(Z — Z 0 )/c*] sin# (69) 




Zo: Axial Distance to Center of Rotations 


(b) Journal rotational motions 


Figure 7. Description of journal motions 

in an annular seal 
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RA8„ 


kl 6 ^ A ^ COs9 + As r sin °- :i - : ^ L (Z-Zo)sin0+ * 


RA6„ 


and 


(Z - Zq ) cos 9] (70) 

Zo = z„/r, T = uti i = x /rf (n) 

where 6 = z = x/R is the dimensionless circumferentia| coordinal( , - n (he 
unwrapped fluid film plane (refer to Figs. 1 to 3), Z = Z/R the dimensionless 
inertial axial coordinate, and *(,) is a general function that describes the 
clearance variations in the axial coordinate. 

For small amplitude motions, all the dependent variables as well 

- the fluid properties (pMetc.) are expressed as the superposition of zero.h- 

and first-order fields representing the steady state and dynamic motion conditions, 
respectively. In general, we let 

<f> = <f>o + e iT (A^ r +Ae<& + RA 8 X, flAA . 

™ ^ + + = *« + *' r (72) 

where 

< t > -u,v,p,f,p,fi,/3 u l ex ,ky,...,etc. and j = X, Y, 6 X , S Y (73) 

Expansion of the governing equations in the perturbation variables yields the 
dimensionless zeroth- and first-order flow equations presented in Appendix C. 

Boundary Conditions for the Flow 

The boundary conditions for both the zeroth-order and the first-order flow 
variables are given and discussed for the three different configurations. The first 
set refers to a circular (360°) bearing, the second one to a characteristic bearing 
pad, and the third set to an annular pressure seal. 

B oundary Conditions for a circular f36Q°'| bearing 

(a) On the 360-degree extended film land, the pressure, velocity, and tern- 
perature fields are continuous and single-valued in the circumferential direction: 

“o,fo ,po,fo(i,S) = u 0 ,»o,po,fo(i -I- 2ir, y) 


(74) 



48 


u h v jiPi>Tj(x,y) - Uj,Vj,pj,Tj(x + 2n,y) (75) 

(b) At the bearing side discharge planes, the fluid pressure is equal to 

specified values of discharge or sump pressures if no end restrictions are present, 
i.e., 

Po(x,Ir) = p R (x) (76) 

p 0 (x,-l L ) = p L (x) (77) 

In general the discharge pressures are uniform and constant. However, in 
some cryogenic turbopump applications the bearing may be located close to the 
pump-impeller discharge. In this case, the sump pressures are non-uniform but 

rotationally symmetric. The model expresses the discharge pressures in a series 
form given as: 

Pr = Pro + X] P Rc cos ( n ^) + Pr 3 sin (n$) (78) 

71= 1 

Pl = Plo + XI P Lc cos ( n 0) + p\ 3 (79) 

71=1 

where the p £, p" terms correspond to the Fourier coefficients of a series expansion 
of the non-uniform discharge (sump) pressures. 

A value of dynamic pressure ( Pj , j = X, Y, S X ,S Y ) equal to zero is imposed at 
all locations where time invariant pressures are specified, i.e. 

Pj(x,l R ) = 0 (80) 

Pj(x,-l L ) = 0. (81) 

On the other hand, if end seals or wear rings are present at the discharge 
sides of the bearings and the fluid flows out of the bearing to the sump conditions 
[V • n > 0), then there is a pressure increase (ram pressure) at the discharge sides: 

Po(x,Ir) = P R (x) + - Po IC R R e *vZ\ lR ' i , if t’oliR.i > 0 


(82) 
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Po {x,-l L ) = p L (x) + -p 0 K. L R e * p vl\- lL>i ., if y 0 |_, t)£ < 0 (83) 

Pj (x,l R ) = ^/CrR^vo^PoVj +p jVo )\ lRii , if t , 0 | iR i > o ( 84 ) 

P](*,-h) = ^LRe p v 0 (2p 0 v j + pjVo)\-i Li£ , if v 0 |-/ t ,i < 0 (85) 

where K r ,Kl are the e nd seal discharge coefficients at the right and left sides 
of the bearing, respectively. Typical values for K range from 1.0 to 4.0 (Childs, 
1989). However, applications such as those with wear end rings of tight clearances 
present larger values of K's (Scharrer et al., 1992b). 

The temperatures at the discharge planes can be obtained from their upwind 
values on the film lands (Roach, 1976). 

Boundary C onditions for a Bearing Pad 

At low rotational speeds, the pressures at the side boundaries of the pad 
(V = ~l L & Ir) are essentially constant and equal to specified values of ambient 
or sump pressures, similar to those for a circular (360°) bearing. A linear 
pressure distribution is assumed along the leading and trailing edges of the pads 
(downstream and upstream edges of the grooves). 

On the other hand, significant momentum changes occur at the pad leading 
edge at high journal surface speeds. The fluid entering the filmlands at the leading 
edge is assumed to have a dynamic head equal to some fraction of a reference 
dynamic pressure based on the bearing surface speed (Burton and Carper, 1967; 
Smalley et al., 1974; Mori et al., 1991; Ettles and Cameron, 1968), i.e., 

&Po{xinlet,y) = -/CpP°i? e * ( — ) (86) 


Pj(xinlet,y) = 0 


(87) 



The coefficient (K P ) ; s fln . . 50 

Burton and Carper ( 1967 ) su ? (eXPenmental ) ram Pressure factor. 

With large turbulent levels. ” & ^ ^ ^ = ^ ** Wgh Speed fl ™s 

Across the line, ~l L < o < / ,, 

_ „ y < l R, the pad leading edee th P fl 

as P^llel to the bearing Dad : 8 g k fl ° W 15 re garded 

8 Pad ’ ,<e -’ the ax ^ a l velocity is null: 


V °( X inlet, y) = 0 
V j( x inlet,y ) = 0 

Temperatures at the leading and trailing edges of th n 
global -e. , 6 ages of the pads 


en ergy balance equation 


before. 


at a 


( 88 ) 


(89) 


are given by the 


groove for different flow conditions 


as shown 


for an annular pressure J , ^ ^ 
-al inlet and exit p , a „es f ollow: ( ”' B ° Undary - ‘he 

(a) Due to fluid inertia, the local a l 

stagnant conditions at the upstream seal ^7 ^ 7^1 ^ 

-let causes a sudden pressure drop. The entra ^ '' ° C " y “ ‘ he ^ 

(5=0) is modeled b y a simn, R P """ * th ' ^ ^ plane 

J 3 SlmpIe Bernoulli equation: 

;f fn> _1 

V 2 (1 +t)Re p (90) 

Pn =?,(*,<>) = - aMiVji+v ^ . foM>o 

The rnlet circumferential velocities are given by 


U 01 = u 0 (x,0) = aA; 


u n = «j(*,o) = o 


( 92 ) 
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where a is a pre-swirl factor, and “1” denotes the seal entrance location. 

The entrance bulk-flow temperatures are equal to (Yang et al., 1993a) 

fn = fo(i ’ 0) = 1 - f(‘+{)Nl-Afcl£„ (93) 

Tn = + if „ ol>0 _ (94) 

For incompressible liquids, jj, = 0 and p n = 1 , therefore, 


f oi = 1 and Tji = 0. ( 95 ) 

For ideal gases, 1 - & = 0, therefore, 

fo,=f.(x,0) = l-^ (1 + fR and f n = ~ il+()EcV0lVjl „ t , oi>0 

(96) 

(b) At the seal exit plane (y = L/R), the discharge pressure for subsonic 
conditions is: 

Po(x,L/R)=p r ; P j(x,L/R) = o (97) 

If an end seal or wear ring is present, the discharge pressure can be treated in the 
same way as for end seals in a HJB. 

Other Boundary Conditions 


If the bearing is symmetric both in geometry and operating conditions, then, 
the axial velocity and the axial gradients of all the primitive variables are zero at 
the circumferential central plane (y = 0) of the bearing: 


<2 

© 

'Hi 

O 

II 

o 

(98) 

o 

II 

o 
1 H 

(99) 


duo dpo dT ( ) _ 


(100) 


Heat Transfer Between Fluid Film and Bounding Solids 

p rature boundary conditions at the fluid-journal and fluid-bushing 
interfaces are modeled by heat transfer to the journal and the bearing (Q. by 

Eq.(12)) and directly incorporated into the energy equation (1!) (Constantinescu 
1973). 


An adiabatic flow process ( 0 .- 0 ) is shown to be a reasonable approximation 
for turbulent flow HJBs and annular seals with large pressure drops across 
eal/beanng length. In such conditions, the heat generated in the fluid 
is carried away mainly by fluid advection while heat conduction through the 
g/ tator is relatively small. The adiabatic flow assumption greatly simplifies 

the analysis, and numerical predictions based on this assumption correlate well 
with experimental data (Yang et al., 1993a,b). 

Regarding pad journal bearings without external pressurization, the adiabatic 

flow assumption is highly conservative. Without significant supply flow of fresh 

fluid into the fluid film, hot lubricant recirculates in the journal bearing and 

causes a much larger temperature increase than in a HJB or an annular pressure 

seal. Heat convection by the side (axial) flow is greatly reduced and heat transfer 

between the fluid film and the bounding solids is significant for this bearing 
configuration. 


Different THD models exist to address the heat transfer from the fluid film 
to the bounding solids. In a full THD analysis, the solution of the Reynolds 
equation (or the continuity and the Navier-Stokes equations) is coupled to the 
energy equation in the fluid film and the heat conduction equations in the solids. 
This coupling leads to a nonlinear iterative problem which is costly and may be 
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sensitive and prone to numerical instabilities. Another shortcoming of the full 
THD treatment is that the bearing pad is never an isolated, single component 
but part of an assembly that includes bearing liner, housing, pedestals, and seals, 
often made of different materials (Pinkus, 1990). Therefore, a full THD solution 

is strongly machine dependent and not considered necessary in the context of the 
current analysis. 

More simplified heat conduction models have been used to treat the heat flow 
in the journal and bushing. Safar and Szeri (1974) considered the journal to be 
isothermal at the mean temperature of the bushing inner surface. Radial heat flow 
into the bushing of a hydrodynamic journal bearing was considered as a reasonable 
approximation ut this might lead to error in the calculated maximum temperature 
in the fluid film. This over-simplification allows calculation of the global heat flow 
from the heat source (fluid film) to the heat sink (ambient) even for composite 
bearing assemblies. Suganami and Szeri (1979) found that at high rotational 
speeds an adiabatic journal surface assumption might be more appropriate than 
the isothermal surface since heat generation in the fluid film is large in comparison 

with the heat capacity of the journal and the time available to absorb this energy 
is short. 

In this analysis and the computer code being developed, different simplified 
THD models are adapted (Fig. 8). Options included are: 

(a) Adiabatic journal and bearing surfaces Q 3 = 0, 

(b) Adiabatic bearing and isothermal journal surfaces, Q, = hj(T — Tj ), 

(c) Adiabatic journal and isothermal bearing surfaces, Q s = hs(T — Tb), and 

(d) Isothermal journal and bearing surfaces (T = T, uppiy ). 

A more detailed model including radial heat transfer through the sleeve or 

stator will be a subject of investigation on Phase II (1994) of the present work. 
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(c) Adiabatic Journal & Iso- 
thermal Bearing Surfaces 


(d) Isothermal Journal 
& Bearing Surfaces 


Fig, 8 Thermal Models of Heat Transfer with 
Different Surface Boundary Conditions 
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Film Forces, Torque, and Dynamic Force Coefficients 


Once the zeroth-order and the first-order equations are solved, the static and 

dynamic performance characteristics of the bearing can be evaluated. Fluid film 

forces and moments are calculated by integration of the pressure field over the pad 

surfaces. The components of the static equilibrium force and moment are given 
by: 


and 


Fxo 
Fy o 


= APR 2 



cos 6 
sin# 


dOdy 


( 102 ) 


Mxo 

Myo 


= APR 3 



p !(2 - Zo) 


— sin 6 
cos 9 


dOdy 


( 103 ) 


where corresponds to the zeroth-order pressure field for the k-th pad. 

The perturbation analysis allows the dynamic force and moment coefficients 
due to journal center displacements and journal axis rotations to be obtained from 
the general expression for dynamic force given as: 
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The dynamic-force coefficients defined by E,.(104) are important measures of 
dynamic bearing performance since they influence the critical speeds, the resonant 
amplitude response, and rotordynamic stability of the rotor-bearing system. 
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The above expression (Eq.(104)) allows the dynamic force coefficients to be 
calculated from integration of the first-order complex pressure field ( Pi ) over the 
pad surfaces. The force coefficients due to journal center displacements are defined 
by the following relationships (San Andres, 1993a): 


Kij a; M-ij + = — A.P 


R 2 ^ f l * , 

~ \ / pMdddy ; ij =X,Y 

* k- rl J-IlJO J 


(105) 


The force coefficients due to journal axis rotations are equal to 

-*** £ /" r Pi *< « = * r (IM , 

* J ~l rJ 0 3 


The moment coefficients about Z 0 due to journal center displacements 
given by: 


are 


r , 2 , a- . i? 3 ‘^l± d r l * r $ pad 

6,i ~ W = -&P— Elf pHZ-Zo)^, d0dy; i,j = X,Y 

* k=l J -ll J o J 


(107) 


and, the moment coefficients due to rotor axis rotations about Z„ are equal 


to 


RsiS, w Ms t = —A P 


R*_^ d r l * r e 

c * k=1 J-iJo 


Pg.(Z—Zo)hs i dddy\ 




where 


i,j = X,Y 
(108) 


— cos 9, hy = —hs x = sin# (109) 

It is noted that the characterization of force coefficients in terms of stiffness, 

damping and inertia parameters is appropriate for incompressible fluids only. 

For highly compressible liquids and gases, the force coefficients are themselves 
functions of frequency. 
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(f) stability indicators such as whirl frequency ratios for lateral and conical 
journal motions, and equivalent bearing stiffness, 

,g) complete pressure, velocity, and temperature fields on the flow region, as well 
as density and viscosity distributions. 


The hydros ealt has been developed in standard Fortran 77 in a SGI 4D- 
35 computer. A detailed User’s Manual and Tutorial accompany the hydrosealt 
program (San Andres, 1993d) A great number of cases has been run to test 
the hydrosealt code. Comparisons show that for L0 2 damper seals, results 
from hydrosealt code are the same as those from a especialized seal code (Yang 
et al., 1993a). For the single-row, five-recess water HJBs tested at Texas A&M 
ity, predictions from the hydrosealt code show good agreement with the 


results from experiments. Refer to the Examples Manual of hydrosealt for a 
description and discussion of the examples presented (San Andres, 1993e). 
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RESULTS AND DISCUSSION 

The analysis presents a general thermal flow model for the static and dynamic 
force and moment response in fully developed, variable properties, turbulent flow 
fluid film bearings. The bearing geometries of interest correspond to hydrostatic 

journal pad bearings (HJBs), annular pressure seals or damper bearing seals, and 
fixed arc pad journal bearings. 

Numerical predictions from the present THD model and prior analytical de- 
velopments have been presented in several earlier publications. The analytical ef- 
fort has always stressed in presenting comparisons with experimental results when 
these have been available in the open literature. San Andres et al.( 1991c, 1993b) 
and Yang et al.(1993a,b) present results for the effects of journal misalignment and 
thermal effects m annular pressure seals, respectively. Kurtin et al.(1993), Mosher 
(1993) and Franchek et al.(1993) present a direct comparison between experimen- 
tal measurements of the static and dynamic force response of a 5 recess, water 
hydrostatic bearing with numerical predictions from the present model. Yang et 
al. (1993c, d) also present comparisons between the THD adiabatic and isothermal 
models predictions with experimental results for turbulent flow liquid hydrogen 
and water HJBs, respectively. In all cases the correlations with the experimental 
measurements are reported to be good and provide validity for the analysis. 

This section presents results for several relevant fluid film bearing applica- 
tions. Many more examples are discussed in the hydrosealt Examples Manual (San 
Andres, 1993e). The first example studied refers to a direct comparison between 
the load capacity and rotordynamic force coefficients for a hydrostatic bearing 
and a damper bearing seal operating in liquid oxygen. The HJB and damper seal 
have been designed to replace the duplex ball bearings next to the left inducer in 
the liquid oxygen high pressure turbopump ( L0 X HPOTP). The fluid operating 
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conditions (pressure and temperature) as well as the actual bearing clearance, and 
most importantly, the load supported by the bearings are a function of the rotating 
speed of the pump. Information relevant to the load characteristics in the HOPTP 
were obtained from Shoup (1993) and the fluid operating conditions directly ex- 
tracted from the technical report of Heshmat (1991). The hydrostatic bearing 
has 6 recesses of rectangular shape and orifice restrictors, while the damper seal 
consists of two parallel annular seals of convergent tapered clearance and sepa- 
rated by a deep feeding central groove (Von Prageneu, 1990). The damper seal 
has a rough stator surface of the knurled type while the hydrostatic bearing and 
journal surfaces are perfectly smooth. This type of damper seal geometry is also 
known as an annular hydrostatic bearing (Scharrer et al., 1992c). Table 1 shows 
a description of the bearing geometries, and it also presents the actual clearances, 
supply and discharge pressures and supply temperature for liquid oxygen, and 
he nominal load acting on the bearings as the operating speed increases from 
14,035 cpm to 30,367 cpm. Note that the load and pressures are proportional to 
the second power of the rotational speed. At the nominal operating conditions, 
here taken as 26,000 cpm, the nominal clearance in the HJB is equal to 0.175 
mm, while the inlet and exit clearances in the damper seal are equal to 0.221 mm 

d 0.129 mm (ratio — 1.715) with an average clearance identical to that of the 
hydrostatic journal bearing. 

The hydrostatic bearing is designed for operation at the nominal speed with 
a concentric pressure ratio equal to 0.60 to provide maximum direct stiffness 
coefficients. On the other hand, the ratio of inlet to discharge clearance in 
the damper seal has been optimized to also obtain the largest direct stiffness 
coefficients. The maximum specific load (load divided by bearing projected area) 
is equal to 6.55 MPa (950 psi) and 7.22 MPa (1048 psi) for the hydrostatic bearing 


Table 1 . Geometry and operating characteristics of hydrostatic bearing 
and parallel damper seals for load support in LOx HPOTP. 

Hydrostatic Bearing- Number of recesses Nrec = 6 

Diameter fm 
85.1 mm 
(3.35 in) 


Length £L) 
48.85 mm 
(1.92 in) 


Clearance (H*) 
1 75.2pm 
(6.9 mils) 


Axial L engthf/) Recess Are 
24.42 mm 30 deg 
(0.96 in) 


Recess Denth 
508pm 
(20 mils) 

journal and bearing surface conditions: smooth. 

UD=0. 57 ; Hr/C =2 . 9i //L=o.5, Recess area ratio=0.25, C/R =0.0041 

Orifice Cd=0.90; diameter do = 4.44 mm 
Recess edge coefficients £xu=0.0; £xd=0.50; ^y=0.0 

PamDer Seal. Two parallel seal land (2 x L) (groove width unspecified) 

~851 t mm D) ^22 2^ Clearance Rati o Clearance 

22 2 mm 221 - 3 nm 129.1pm 175.2pm 

(3.35 in) (0.874 in) (8.7 mils) (5.08 mils) (6.9 mils) 

journal smooth, bearing rough, r/h = 0.044 (Knurled) 

Entrance coefficient ^y=0.25 

Operating Condition*; 


Speed 

Psupply 

Pa 

Tsupply 

Tsat 

Load 

N 

com 

MPa 


K 

K 

14,035 

16.00 

1.654 

102.77 

128.9 

5,828 

19,732 

29.59 

1.792 

107.22 

130.5 

11,519 

26,000 

39.60 

2.089 

110.55 

133.7 

20,000 

28,340 

43.60 

2.434 

111.33 

137.0 

23,762 

30,367 

55.69 

2.551 

115.00 

138.0 

27,282 


Reynolds# 
Rec 
72,816 
100,700 
130,130 
143,720 
147,600 


Average 
Clearance* 
im) 

180 
177.7 

175.2 

174.3 
173.2 


Inlet swirl ratio cx= 0.50 

Fluid: L02 (liquid oxygen) at 1 1 0.6 K (200R) 
-£- ( MPa > 6_(kg/m3) u (E-3 Pa s) 


39.60 

2.09 


1134.4 

1038.1 


0.177 

0.124 


1,606.4 

1,827.2 


fhydrostatic 


Jam per seal 
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and damper seal, respectively. These specific loads are very large considering the 
nature of the bearing application with a fluid of very low viscosity such as liquid 
oxygen. The calculations for the static and dynamic performance characterises 
of the HJB and seal are performed using the adiabatic journal (rotor) and bearing 
(stator) thermal model. Numerical calculations for the damper seal are performed 
only for the thin land portion and then multiplied by two. The central groove and 
its effect on load capacity and dynamic force coefficients are altogether neglected. 
This oversimplification seems appropriate as a first attempt to correlate the 
performance of the HJB and seal. However, it is now known that central feeding 
grooves do have a pronounced effect on the dynamic force response of fluid film 
bearings, in particular with regard to inertia and cross-coupled damping force 
coefficients (Arauz et al., 1993, Lindsey, 1993). 

Figure 9 depicts the fluid supply temperature, bearing load and pressure 
drop across the HJB and seals as the rotational speed of the pump increases. The 
largest load of 27,282 N (6,137 lbs) corresponds to the highest operating speed. 
The Figure also shows the values of the nominal circumfererential flow Reynolds 
number Re c =p^Rflc m / /i* to range from approximately 70,000 to 150,000 as the 
rotor speed rises. The flow in the bearings is then turbulent. Figure 10 presents 
the HJB and damper seal flow rate and drag torque as the journal speed increases. 
The hydrostatic bearing (h) shows approximately 14 percent more flow rate than 
the damper seal (s), while it produces approximately 27 percent less drag torque 
than the seal at the largest operating speed. These results are a direct consequence 
of the rough stator surface in the damper seal. 

Figure 11 presents the journal operating eccentricity and attitude angle as 
well as the maximum temperature rise in the HJB and damper seal as the 
operating speed increases. Note that as the speed rises so does the applied load 



1 40 I- LOX: 6 recess HJB and damper seal 
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Fig. 9. Supply temperature, pressure drop, Reynolds number 

and load vs. rotational speed for LOx 6 recess HJB and damper seal. 




6recess HJB and damper seal 
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Fig.10. Flow rate and drag torque vs. rotational speed for 
LOx (h) 6 recess HJB and (s) parallel damper seal. 
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at a rate proportional to the journal speed squared. The load for the HJB is 
directed towards the bottom of a bearing recess (X direction). The dimensionless 
journal eccentricity has been determined as the ratio between the journal offcenter 
displacement divided by the nominal clearance at 26,000 cpm, i.e. 0. 1 75 mm. Note 
that the attitude angle is less than 10° for both bearings and the journal eccentric 
displacement is rather moderate considering the magnitude of the loads applied. 
The low value in the attitude angle indicates dominance of hydrostatic effects 
over hydrodynamic effects. The maximum temperature in the film lands of the 
hydrostatic bearing and seals increases rapidly with journal speed. The HJB shows 
a larger thermal differential between the bearing supply condition and discharge 
planes. The results also indicate that both HJB and damper seal operate well 
below the critical LO x temperature of T c = 154.6 °K (278 °R) 

Figures 12 to 14 present the synchronous stiffness, damping and inertia force 
coefficients versus the rotational speed for the HJB and damper seals, respectively. 
Note the similar values between all direct coefficients (say K xx and K „ ) and cross 
coupled coefficients (say K„ and -K rx ) which denote that the HJB and damper 
bearing seal have very uniform force coefficients as the operating eccentricity (and 
load) increases with the operating speed. Figure 12 shows a HJB with slightly 
smaller direct stiffness coefficients (A'„ and K„) than the damper seal, while its 
ross coupled stiffness ih xr ) is larger. Note the dominance of hydrostatic (direct) 
coefficients over the cross-coupled coefficients induced solely by journal rotation. 

Figure 13 shows the HJB to have larger direct damping coefficients (C„ and 
C„) than the damper seals (approximately 47 percent higher). These results 
produce at the largest journal speed a whirl frequency ratio (WFR) equal to 0.39 
and 0.45 for the HJB and damper seal, respectively. Thus, in this example when 
the inlet swirl ratio is equal to 0.50, the HJB offers slightly better dynamic stability 


(h): 6 recess HJB and (s): damper seal, fluid:LOX 



Q|6uv ‘si-xeiuj^ ‘Ajiojjjueooe |eujno[ 


Fig.11. Journal eccentricity and attitude angle and maximum temperature 
rise vs. rotational speed for 

LOx (h) 6 recess HJB and (s) parallel damper seal. 



6 recess LOX HJB and damper seal 
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F'9-12. Stiffness Coefficients vs. rotational speed for 

LOx (h) 6 recess HJB and (s) parallel damper seal. 
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characteristics than the parallel lands damper seal. 

Figure 14 shows the damper seal direct inertia force coefficients (M xx and 
M„) to be the on the order of 1.5 kg and large in relationship to the HJB direct 
inertia coefficients. The reason for this lies on the tapered geometry of the damper 
bearing seal with a small clearance at the exit plane. The cross-coupled inertia 
coefficients (Af„ and M rx ) are small in nature. Here it is noted that in practice 
the direct inertia coefficients for the damper seal will be much larger than the ones 
predicted. The central feeding groove will act as a parallel inertia and compliance 
to the seal lands, it will increase the value of the inertia coefficients, and reduce 
substantially the direct dynamic stiffness of the damper seal. Experimental results 
showing this behavior have been given recently by Lindsey (1993). 

The example demonstrates that a properly designed HJB or a parallel damper 
bearing seal can support easily the loads expected in the LO, HP turbopump. 
The results show that the HJB and seal will operate at low eccentricity ratios 
With uniform force coefficients. The case studied shows an important application 
of the hydrostatic principle where good engineering practice utilizes the available 
P ure differential in a cryogenic turbopump to provide a reliable fluid film 
bearing support. It i, worth noting that the original test case included a three 


pad journal bearing with a clearance similar to that of the HJB but with a 
bearing preload of 0.076 mm. The results are not reproduced here because the 
j urnal bearing offered a very low load capacity in comparison with the HJB and 
damper seal, and it also produced a rather large temperature rise. Details of these 
calculations can be found elsewhere (San Andres, 1993e). 


Adams et al. (1992) presented experimental results for the force coefficients 
of a four pad, one recess/pad, hydrostatic bearing operating at low speeds and 
using motor oil SAE 30 at (T = 38 V) as the working fluid. Table 2 shows the 



LOX: 6 recess HJB and damper seal 
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Fig.13. Damping Coefficients vs. rotational speed for 

LOx (h) 6 recess HJB and (s) parallel damper seal. 
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geometry and operating characteristics of the HJB test article with four recesses 
of large area and depth. The maximum operating speed is equal to 2,000 cpm and 
the circumferential flow Reynolds number is equal to 23.0. The flow conditions 
clearly specify a laminar flow regime without significant fluid inertia effects at 
film lands and recess edge on the forced response of the bearing. Adams et al.’s 
tests derive from an experimental research program whose results may be used 
for extrapolation of the dynamic force performance in high-speed, turbulent flow 
hydrostatic bearings in cryogenic turbopumps. Adams et al. present the force 
coefficients in graphical form versus the amplitude of a dynamic eccentricity 
used in the experimental procedure. The numerical predictions are based on 

the adiabatic thermal model with the journal center operating at the concentric 
position. 

Table 3 presents a comparison between the predicted values and the experi- 
mental force coefficients determined by Adams et al. The values for the measure- 
ments correspond to average values along with the largest positive and negative 
deviations from the average number. The numerical predictions show a good 
agreement with the experimental results although the latter show very wide vari- 
ations. Very similar numerical results were obtained for a flow model with isother- 
mal characteristics. The hydrostatic bearing of Adams et al. (1992) demonstrates 
that the bearing flow rate and direct stiffness coefficients are not affected by ro- 
tational speed. On the other hand, the cross-coupled stiffness coefficients are 
directly proportional to journal speed. The results are representative of laminar 
flow hydrostatic bearings and can not be extended to turbulent flow operating 
conditions. Inertia coefficients are not reported here. Added mass coefficients are 
difficult to calculate (and measure) accurately for small values of the squeeze film 
Reynolds number (Re,). The experimental results for inertia force coefficients 



Tab ' e 2 ‘ °!°r'T and ° pera,in 9 Conditions of Four Pad Hydrostatic 
Bearing from Adams et al. (1992). 

Number of Pads: Npad = 4 , Number recesses/pad Nrec = 1 


Pads: circumferential length ©p=80.50 deg, groove width = 9.50deg. 

Recess: circumferential length 0r=56.deg 

1^0=0.47, f/L=0.659, 0r/@p=o.696 , Area Ratio=0.459, C/R=0.0037, Hr/C=22.6 
journal and bearing surface conditions: smooth. 

Orifice Cd=1 .00, diameter d,=,. 592mm , or concentric recess presure ratio pr=0.40, 
Recess edge coefficients £xu=0.0; £xd=0.00; £ y =0.0, 

Pad leading edge coefficient, k=0.64 

Operatin g Parameters- 

Rotational speed: 1,000 and 2,000 cpm 

Pressure supply, Ps=1 ,206 MPa (175psig), 2.586 MPa (375psig) 
exit, Pa= 0.00 MPa ( Opsig) ' p 9) 

Recess Pressure, Pr=0.483 MPa (70psig), 1.034 MPa (ISOpsig) 

Fluid: SAE 30oil at 310K (100F),ps=0.09787Pa-s, P s=p a =890.00 Kg/m3 

TYP Reyn olds numbers • 

Redrc=pnRc/p s= i 1 .5 (Ikcpm) , 23.0 (at 2kcpm) 
eaxial m/2ttDR pa)=4.03 ( Pr=70psig), 8.64 (Pr.150psig) LAMINAR FLOW 
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seem to be erroneous with identified magnitudes as large as 100 kg. Figure 15 
Shows the calculated pressure field a. the concentric position for the four pad 
hydrostatic bearing studied. The figure depicts a pressure field typical of a pure 
hydrostatic bearing without fluid inertia effects. 


y ynamic journal bearings, in particular those of a compliant bearing 

surface nature, are currently a subject of attention for cryogenic applications 

(Heshmat, 1991). Experimental results for static and dynamic force performance 

characteristics in pure hydrodynamic journal bearings with cryogenic liquids have 

eported in the open literature. The interested reader needs to refer 

to the literature review provided by San Andres (1993f) on the subject of foil 
bearings. 


The present THD model for journal bearings is used to analyse an axially 
grooved journal bearing as tested by Tonnesen and Hansen (1981). The published 
experimental measurements are very complete and detailed. This test bearing has 
then become an ideal example to validate THD models by several researchers (see 
Knight and Niewiarowski, 1990). Details of the bearing geometry and operating 
conditions are given in Table 4. The calculations are performed at a journal 
speed of 3,600 cpm with external loads directed towards the bottom pad and 
equal to 5,600 N and 8,600 N. The bearing analyzed handles a viscous lubricant 
and the smallness of the film clearance and rotational speed determine laminar 
flow conditions. For the calculations the adiabatic flow assumption was used 
although this is too crude since the measurements definitely show there is heat 
flow through the bearing surfaces and rotating journal. A groove mixing coefficient 
(A) equal to 0.80 was used on the calculations as derived from the experimental 
results. Numerical predictions for the bearing centerline dimensionless pressure 
and temperature are shown in Figures 16 and 17, respectively. Table 4 also 



Table 3. Performance Characteristics of Four Pad Hydrostatic Bearing at 

eccentricity ex=0.0. Comparison to test results from Adams et al.(1992) 


Torque temperature 


Reynolds # 
Min/Max 
2.66/9.59 


Squeeze 
Reynolds Re s 
0.042 


Speed 

KcDm 

Precess 

Mpa 

KXX 

Kyy 

MN/m 

KXY 

kyx 


1.0 

0.48 

26.28 

21.17+2.46 

-1.97 

26.28 

22.32+3.28 
-1 .66 

4.06 

5.25+1.15 

-0.82 

-4.07 

-3.05+1.28 

-0.88 

predicted 

experimental 


1.03 

56.16 

44.81+8.86 

-5.42 

56.20 

43.66+5.58 

-3.28 

4.94 

5.91+1.97 

-2.95 

-4.78 

-4.27+3.28 
-2 62 

predicted 

experimental 

2.0 

0.48 

26.00 

19.17+1.84 

-1.84 

26.00 

22.76+2.45 

-2.80 

8.19 

7.88+1.57 

-1.58 

-8.19 

-7.35+1.57 

-2.10 

predicted 

experimental 


1.03 

55.76 

A A tzr\ . O a 

55.76 

7.99 

-7.98 

predicted 


9.80+3.85 -9.80+4.55 experimental 

^-3.50 -4.90 


Damping Coefficients 


Speed 
Kcpm 
1.0 


Precess 
Mpa 
0.48 


cxx cyy 

KNs/m 
78.08 79.07 

65.65+10.10 70.70+10.00 


CXY CYX 

-0.29 predicted 

-4.33+4.32 experimental 

-8.65 

5.98 predicted 

-30.25+20.80 experimental 
-82.30 

-0.40 predicted 

-0.21+6.29 experimental 
^ 3.43 

predicted 
experimental 


Experimental values obtained from graphical values shown by Adams et al (1992) 
Values refer to average quantities with maximum and minimum deviations from average. 


NPadS=4 1 rCCeSS/pad ’ ‘■-IH.nrnm. 0=53. 97mm, l=35.56mm, 

groove width=9 525 mm R ' S6 ' 134mm ’ c=210.8mm, Hr/c=22.S, 

v ve mam 9.525mm, recess Area ratio=0.4590 


P S 


90 


180 


Recess Pressure Ratio- 


270 


0 


ill 


groove 


\P/Ps 


Recess 


land 


Ps-1'206MPa(17Spsig), PR=0.482MPa(70psig), 1000 rpm 
Qmd: on SAE 30 , 37.5C (100F), ,-0.097Pa s, ,=890kg,m3 

Rec=11.50, Re™, =13.37, 0=0.286 kg/s, do=1.59mm, Cd=1.0 

Figure 15. Pressure Field m a 4 recess Hydrostatic Pad Bearing 
oncentric Position. Adams et ai. (1992) bearing 
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Table 4. Geometry and operating characteristics of axially grooved 
journal bearing from Tonnesen and Hansen (1981). 

2 pad journal bearing:. 

Le " q ' hfL) C 'lnarance (C*) Part a rc Grnnvn width 

55 "™ 69.0 pm 17tTdeg iodSS 

(3.93 in) (2.16 in) (2.71 mils) 

journal and bearing surface conditions: smooth,: L7D=0.55: C/R=0.0014 

Thermal mixing groove coefficient X= 0.80 

011 “ T! • e i mpera,ure 50 C (122 F) and Psupply=0.1 4MPa (20.3psig) 
P (kq/m3) li l Pa.sl Cnll/kn-K' ’ P 9) 

50 850 0.0182 2,000 

90 850 0.0056 2,000 

Operating Condition: journal speed 3,500 cpm 


Load = 5,600 N 

Case e/c 

THD soln. 0.53 


Min T (C) Max T (Cl Max P (MPat 


THD soln. 0.53 63.22 64.50 

: 55.60 7223 

experimental 0.56 58.00 62.00 

57.50 72.00 


0.14 

2.58 


0.14 

2.60 


Pad 

# 1 cavitated 

# 2 loaded 


Load = 8,600 N 

^ — Min T (C) Ma x T (Cl Max P tMPa\ 
THD soln. 0.63 63.10 64.70 0.14 

55.50 75.10 4 31 

experimental 0.66 58.00 noTgiven 0A4 

57.50 74.00 in 


0.14 

4,31 

0.14 

5.10 




Pad 

# 1 cavitated 

# 2 loaded , 






Tonnesen bearing 3.5Krpm; Pa=0.14MPa, Ts=273.3+50C 



s±/0 ^x(si-i) sjn^jadiLiai eujiJajuaQ 


Fig.17. Centerline temperature for axially grooved journal bearing. 
Bearing parameters as per Tonnesen et al. (1981). Loads 
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regards the fluid flow as laminar, although large journal rotational speeds were 
considered in the application studied. The bearing geometry and operating 
conditions for the bearing studied by Braun et al. are given in Table 5. The 
journal bearing has a very small clearance equal to 0.0254 mm (0.001 in). At the 
rotational speed of 14,000 cpm the circumferential flow Reynolds number is Re c 
= 12,027 with fluid properties evaluated at a supply pressure and temperature 
equal to 1.72 Mpa (250 psi) and 111 0 # (200°R), respectively. 

The results presented by Braun et al. are most peculiar since they show 
a temperature field (Fig. 5 of cited reference) which is a maximum at the 
bearing center line and decreasing towards the axial discharge planes. Also 
the temperature field does not show a zero temperature gradient at the bearing 
centerplane as theory requires. Another shortcoming for the referred THD analysis 
is related to the adiabatic (journal and bearing) model used for a full 360° journal 
bearing as the calculations suggest. In a continous journal bearing without any 
feeding grooves, the adiabatic model is physically incorrect since, as the flow 
equations also show, the fluid at the bearing circumferential middle plane is 
’’trapped” because the axial pressure gradient is null at this location. For this 
flow condition, a fluid particle at the midplane rotates ad-infinitum around the 
bearing and its temperature should increase without bound. These statements are 
by no means trivial since they were discovered (and thought over) after a long and 
tedious unsuccesful attempt to reproduce the results given in the cited reference. 

For the calculations presented here it was decided to keep the adiabatic model 
(insulated journal and bearing surfaces) but to have the journal bearing with a 
feeding groove of arc length 10° at an angular location equal to 70° as shown on 
Table 5. The thermal mixing groove coefficient (A) used in the calculations is 
equal to zero, i.e. only fresh fluid at the supply temperature enters the bearing 
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■ 1 » vi 1 

111.1 1 723 . — 

(200 R) (250psia) ’- 033 


121.9 


[r6 Pa 


I 1 

1,836.8 


l=K) 


at SUpp,y ( ,nlet ) conditions 


. . ^ (iniet) condj 

*-Oa critical temperature and pressure: 154.8 K (2 78 6fl . 
afiSaa0SL£2 fc journal speed 6 000 ' ^ ' 737 ' 4 P * # 

7 74027 at 14 kcpm 


deeding 

groove 


70 deg 


Load 


X 
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downstream of the feeding groove. The value used is rather conservative but highly 
recommended as it would become evident latter. Preliminary calculations with 
the hydrosealt program showed that the temperature within the fluid film region 
increased rapidly and this caused the saturation pressure of the liquid to drop 
sharply. For example, at a temperature of 111 0 * the liquid saturation pressure is 
equal to 0.582 MPa and much lower than the feeding groove pressure (1.723 MPa). 
However, at a film temperature of 126 °K the saturation pressure is only 1.425 MPa 
and very close to the feeding pressure. It is well known that journal bearings are 
able to produce hydrodynamic pressures above (positive) and below (negative) 
the reference supply value. Liquid flashing (dispersion into a gaseous phase) is 
caused since these ’’negative” pressures can not be lower than the liquid saturation 
pressure. Here we will refer to this phenomena as some sort of cavitation since the 
computer model is unable to distinguish between liquid boiling (two-phase flow) 
and actual liquid cavitation. For the 14,000 cpm example, the cavitation pressure 
is set equal to the supply pressure (1.72 MPa) in the results that follow. On the 
other hand, for 6,000 cpm operation, the cavitation pressure is equal to 1.0 Mpa. 

Table 6 provides a summary of the results calculated with the present analysis. 
Two models are considered (a) adiabatic journal and bearing surfaces, and (b) 
isothermal bearing without thermal effects. The calculations are performed for 
two loads equal to 2.5 KN and 5.0 KN. Figure 18 shows the centerline pressure field 
and temperature for the two loads considered. The pressure and temperature fields 
are symmetric about the bearing midplane with the temperature field varying very 
little in the axial direction. At the largest load considered (5.0 KN), the maximum 
temperature predicted within the fluid film is equal to 127.5 °K, i.e. 16.4°AT above 
the supply value, and the operating journal eccentricity and attitude angle are 
equal to 0.371 and 77.2°, respectively. On the other hand, the results of Braun 
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Table 6. LOx hydrodynamic journal bearing. Comparison of results lor two 

Braim e°aJ. S (198n5* ,IC and ls0lhermal models and “"“"aliens from 

Cases 

B Predictions from Braun et al. (1987b). (*) denotes no values available 
A Present model : adiabatic journal and bearing surfaces with variable fluid properties 
isothermal conditions with variable fluid properties 
artificial setting of laminar flow and no fluid inertia 


Load 

N 

Case e/c <f> Torque Max P 

n N-m MPa 

Max T 

°K 

Reynolds # 

max / min 

2,500 

A 

1 

0.191 90.0 1.80 2.52 

_ 0.213 86.2 2.03 2.55 

126.6 

111.1 

9,050 6,013 

7 290 4 799 

5,000 

5.102 

A 

1 

B 

0.371 77.2 1.84 3.44 

0.381 77.0 2.08 3.48 

0.700 66.1 

127.5 10,613 4,919 
111.1 8,298 3,723 

115.7 * 

Load 

N 

Case 

Kxx Kyx Kyy Kxy 
MN/m 

Cxx 

Cyx Cyy Cxy 
KN/s 

2,500 

A 

1 

50.1 -431.1 98.0 465.8 

_ 25.8 -401.8 109.5 504 7 

604.0 

658.1 

-32.9 625.2 205.6 
5.0 587 2 255 5 

5,000 

5.102 

A 

1 

B 

158.9 -416.4 266.6 634.7 
156.8 -354.0 256.4 674.4 

22.9 -255.1 533.6 52? ft 

777.0 

843.5 

* 

96.9 689.9 372.0 
108.4 598.3 391.0 

Load(N) Case 

Mxx Mvx Mxv Myq 

/ (kn) 


2,500 

A 

1 

58.1 -4.0 59.0 8.6 

_ 59.5 0.9 48.6 4 1 



5,000 

A 

1 

59.0 -2.8 53.9 13.4 

_ 60.3 0.2 40 0 9 a 



Journal SDeed: 6 non mm 



Load 

N 

Case 

e/c <}> Torque Max P 
H N-m MPa 

Max T 
°K 

Reynolds # 

mPY / min 

2,300 

2.292 

A 

B 

0.423 94.2 0.47 2.20 119.0 3,959 1,607 

0.700 90.0 * * i ip n * 

Load 

N 

Case 

Kxx Kyx Kyy Kxy 
MN/m 

Cxx 

KN/s 

Cyx Cyy Cxy 


2,292 


9.3 -194.5 4.55 283.6 
0.0 -128.9 0.00 323.4 
_6.9 -134.9 8.40 329 Q 


710.9 -106.6 878.3 33.4 
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et al. analysis show an operating journal ecentricity equal to 0.7 and an attitude 
angle equal to 66.06 0 with a maximum film temperature equal to 115.8°/f. The 
differences noted between the two THD analyses are easily explained since Braun 
et al’s model considers the flow as laminar without fluid inertia effects. Dynamic 
force coefficients are given in Table 6, while Braun et al. only report values for the 
stiffness coefficients. It is also worth noticing that the results from the adiabatic 
and isothermal models are very similar although the thermal field is quite different 
for both models. A groove mixing coefficient (A) larger than zero will cause the 
temperature rise in the fluid film region to be larger than the ones presented. 

Results for operation at a journal speed equal to 6,000 cpm and a load equal 
to 2,300 N are also shown in Table 6. The maximum temperature predicted by 
the present adiabatic model within the fluid film is equal to 119°f£T, i.e. 7.9 °K 
above the supply value, and the operating journal eccentricity ratio and attitude 
angle are equal to 0.423 and 94.2°, respectively. The bearing shows no cavitation 
and the attitude angle is larger than 90° due to the effect of fluid inertia. On the 
other hand, Braun et al. results present an operating journal ecentricity equal 
to 0.70 and an attitude angle equal to 90° with a maximum film temperature 
equal to 112°.ftT for a load equal to 2,292 N. Again, the difference in results is 
explained by the absence of fluid inertia and turbulence flow effects in the cited 
reference of Braun et til. It is also noted that calculations were also performed for 
a journal bearing with laminar flow conditions by artificially reducing the density 
of the liquid. In this case, the numerical results render a load of 2,600 N at 
an eccentricity ratio equal to 0.70 with stiffness coefficients very similar to those 
reported by Braun et al. (Table 2, 1987b). 

The last example presented refers to a three-pad journal bearing operating 
with liquid oxygen. The geometry and operating characteristics for the rigid 
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Pads bearing are given in Table 7. The data has been taken direct* front 

(1991). The example discusses the bearing load performance at a journal 

speed of 29,830 cpm for two different THD models, namely (a) adiabatic or 

msu ated journal and bearing surfaces, and (b) isothermal bearing surface at the 

supply temperature and adiabatic journal. The bearing clearance is rather small 

0.0381 mm) in order for the bearing to provide a substantial load capacity. 

Note that externally pressurized bearings such as HJBs or damper seals do not 

present this res, notion. The larges, load of 10,000 N determines a specific load in 

the bearing equal to 2.52 MPa (367 psi). The hydrodynamic 3-pad bearing test 

to provide some benchmark data for comparisons with numerical 

tions obtained for a bump-type foil bearing. Details on the static and 

dynamic force performance characteristics of the compliant surface bearing are 
given by San Andres (1993f). 

A thermal mixing groove coefficient (A) equal to zero is taken for all calcula- 
tions. This value ,s very conservative but required to keep the fluid temperature 
-thin tolerable limits for the present application. ,n the numerical computa- 
tions the fluid is not allowed to sustain subambien, pressures, i.e. the mimimum 
pressures are equal to the supply groove pressure of 5.52 MPa (800 psi) This 
assumption is physically correct if the bearing pads are of a compliant nature 

Note that the nominal value of the circumferential flow Reynolds number (Re.) 

IS equal to 32,216. 

Figure 19 shows the bearing equilibrium journal eccentricity ratio and atti- 
tude angle as the bearing load increases from zero to a maximum value of 10,000 
N (2,250 lbs). For small loads the attitude angle is larger than 90” denoting the 
effect of fluid inertia. The adiabatic bearing model shows a slightly smaller op- 
erating eccentricity and larger attitude angle than the isothermal bearing model. 
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Table 7. Liquid Oxygen, three pad hydrodynamic journal bearing. 

Geometry and operating conditions. Data from Heshmat ( 1991 ). 

Number of bearing pads: 3 

Diameter (D) Length fL) Clearance (C*) Pad Arc Groove width 
88.9 mm 44.45 mm 38.1 pm 110 deg 10 deg 

(3.50 in) (1.75 in) (1.50 mils) 

journal and bearing surfaces: smooth.: L/D=0.50; C/R=0.00086 

Thermal mixing groove coefficient X= 0.0 


Fluid: LO2 (liquid oxygen) 

IliK) — P (MPa) 0 (kq/m3) u (E-6 Pa.s) CpfJ/kn-KI 

95 5.52 1,129.6 183.2 1,718 4 

(171 R) (800psia) 


at supply (inlet) conditions 


LO 2 critical temperature and pressure: 154.8 K (278.6 R), 5.083 MPa (737.4 psia) 

Operating Condition: journal speed 29,830 cpm. 

Load varies from 0 to 1 0,000 N (2,245 lbs) 

Nominal circumferential flow Reynolds # (Rec) = 32,616 

Cases : (A) adiabatic journal and adiabatic bearing. 

(I) adiabatic journal and isothermal bearing at Ts 



X 
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The journal eccentricity increases rapidly in a non-linear form with the applied 
load. Figure 20 depicts the centerline dimensionless pressure field and tempera- 
ture rise for a bearing load equal to 7,000 N. Note that the pressure field is almost 
identical for both bearing models while the temperature curves show a different 
character. The results show the third bearing pad to be fully unloaded. Figure 
21 presents the maximum temperature rise for both bearing models. The adia- 
batic model shows a dramatic increase in the bearing exit temperature even for 
the unloaded condition. The maximum temperature rise is over 50 °K (90 °R) at 
the largest load. On the other hand, the isothermal bearing model offers a very 
moderate temperature rise which remains almost invariant for all loads applied. 
In the isothermal bearing case the pad surface acts as a sink for thermal energy 
and prevents the excessive rise of temperature within the film lands (see Figure 
20 ). 

Figure 22 shows the maximum pressure developed on the bearing and the 
drag torque as the bearing load increases. The peak pressure is proportional to 
the applied load, and with the isothermal model showing a larger pressure at 
the largest load. The drag torque is lower for the adiabatic model bearing since 
the larger film temperaturs reduce the liquid viscosity. On the other hand, this 
reduction in fluid viscosity due to temperature also causes an increment in the 
flow Reynolds number and enhances the turbulence flow effects. This may explain 
why the adiabatic bearing model produces a smaller bearing eccentricity than 
the isothermal model. Note that the present THD model assumes the bearing 
clearance not to be affected by thermal phenomena. This seems to be a rather 
poor consideration if the fluid temperature is to rise substantially. 

Figures 23 to 25 present the synchronous stiffness ( Kij ), damping ( Cij ) 
and inertia (M%j) force coefficients versus the applied load for the high speed 
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Fig 20. Centerline dimensionless pressure field and Temperature rise 
for LOX 3-pad bearing. Journal speed 29,830 cpm. 

(A) adiabatic bearing, (I) isothermal bearing 




3 pad LOX bearing at 29830 rpm, Ps=Pa=5.52MPa(800psi) 
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Fig 22. Bearing torque and maximum pressure vs. applied load 
for LOX 3-pad bearing. Journal speed 29,830 cpm. 

(A) adiabatic bearing, (I) isothermal bearing 
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Bearing Load at bottom pad (KN) 

Fig 23. Synchronous stiffness coefficients vs. applied load 

7av° X 3_pad bearin 9- Journal speed 29,830 cpm 
(A) adiabatic bearing, (I) isothermal bearing 


3 pad LOX bearing at 29830 rpm 



[w/s'NX] sjuepjjjeoo 6ujdujeQ snouojipuAQ 


Fig 24. Synchronous dampig coefficients vs. applied load 
for LOX 3-pad bearing. Journal speed 29,830 cpm. 
(A) adiabatic bearing, (I) isothermal bearing 


3 pad LOX bearing at 29830 rpm 
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Fig 25. Inertia force coefficients vs. applied load 

for LOX 3-pad bearing. Journal speed 29,830 cpm. 
(A) adiabatic bearing, (I) isothermal bearing 
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larger equivalent stiffness, its critical mass would be lower since the WFR 
is larger than the adiabatic flow model. 

The LO x bearing examples have shown that a cryogenic hydrodynamic 
j rnal bearing is able to support large loads if the bearing clearance is very small 
and its projected area (L x D) is large. However, these same requirements cause 
thermal effects to be of importance on the static and dynamic force performance 
of the bearing. The actual thermal model to be used for a cryogenic hydrodynamic 

bearing still needs to be determined by considering the thermal paths across the 
bearing and journal surfaces. 

To close this section it is noted that no examples have been presented for 
bearing applications with liquid hydrogen. The ommission has been involuntary 
and several relevant cases can be found in the Examples Manual (San Andres 
1993e) and also in the reference of Yang et al. (1993c). However, it is here 
advanced that thermal effects are rather insignificant for externally pressurized 
LH, bearings since the specific heat of this fluid is much larger than that of LO 
On the other hand, fluid compressibility considerations and a Joule-Thompson 
like effect rule the performance of LH 2 bearings. 



Adiabatic Journal/lsothermal stator 
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Fig 26. Whirl frequency ratio (WFR) vs. applied load 

for LOX 3-pad bearing. Journal speed 29,830 cpm. 
(A) adiabatic bearing, (I) isothermal bearing 




3 pad LOX bearing at 29830 rpm 
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Fig 27. Equivalent stiffness coefficient vs. applied load 

for LOX 3-pad bearing. Journal speed 29,830 cpm. 
(A) adiabatic bearing, (I) isothermal bearing 
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CONCLUSIONS AND RECOMMENDATIONS 

A thermohydrodynamic analysis and computer code for the prediction of the 
static and dynamic force response of fluid film bearings for cryogenic applications 
have been developed. The motion of a cryogenic liquid on the thin film annular 
region of a fluid film bearing is described by a set of mass and momentum 
conservation, and energy transport equations for the primitive turbulent bulk- 
flow variables, and accompanied by thermophysical state equations for evaluation 
of the fluid material properties. Zeroth-order equations describe the fluid flow 
field for a journal equilibrium position, while first-order linear equations govern 
the fluid flow for small amplitude journal center translational motions and journal 
axis conical motions. Solution to the zeroth-order flow field equations provides 
the bearing flow rate, load capacity, restoring moments and torque. Solution 
to the first-order equations determines the 32 rotordynamic force and moment 
coefficients due to journal lateral and angular motions. The analysis includes the 
effects of flow turbulence, fluid inertia, liquid compressibility and thermal energy 
transport on the performance of cryogenic liquid bearings. 

The current technological needs call for reliable and resilient fluid film 
bearing designs to provide maximum operating life with optimum rotordynamic 
characteristics at the lowest cost. The analysis and computer code developed 
constitute practical tools for the performance prediction and design of cryogenic 

liquid hydrostatic journal bearings, annular pressure seals, and cylindrical pad 
bearings. 

The numerical predictions from the program developed correlate favorably 
with experimental results available in the literature. The comparisons performed 
demonstrate the generality of the analysis and validate the computational model, 
and also extend the range of applicability of the 


program to conventional bearings 
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handling viscous lubricants. 

The computer program developed named as hydrosealt calculates: 

1) bearing flowrate (seal leakage), 

2) friction torque and temperature rise, 

3) load capacity (fluid film forces) and restoring moments, 

4) rotordynamic force coefficients due to journal center displacements, 

) rotordynamic force coefficients due to journal axis rotations, 

6) rotordynamic moment coefficients due to journal center displacements, 

7 rotordynamic moment coefficients due to journal axis rotations, 

8) Complete pressure and temperature fields on the bearing surface, as well 

as density and viscosity field variations, with ranges of fluid flow 

Reynolds numbers and Mach numbers. 

for isothermal or adiabatic journal and bearing surface conditions. 

as a function of: 

a) journal (rotor) center eccentricity and journal axis misalignment. 

b) inlet specified circumferential pre-swirl velocity. 

c) general clearance function in the axial direction. 

d) mean surface roughness on bearing and journal film lands. 

The fluid properties (density, viscosity and specific heat) are calculated from 

the 32 term Benedict- Webb-Rubin equations of state as given by NBS Standard 

Reference Database 12 (McCarty, 1986) for the following cryogenic liquids: L0 2 , 

LH 2 , LN 2 , and methane. Other fluids included in the program are water, mineral 
oils and air. 

The analytical research and computational program development need urgent 
support from a complimentary experimental research effort directed to measure 
the dynamic force performance of high speed fluid film bearings with flow 
conditions similar to those found in cryogenic environments. Some empirical flow 
parameters need to be identified from the test data. These coefficients, like for 
example orifice discharge coefficients, recess edge entrance loss factors and shear 
friction coefficients, and thermal mixing groove coefficients are of fundamental 
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importance for the accurate prediction of the flow and force dynamics in high- 
speed, turbulent flow fluid film bearings. 

Further research is needed to continue the development of advanced analysis 
and computer codes for the prediction of the force response of other fluid 
film bearing geometries for cryogenic applications. The study should consider 
extension of the analytical models and program developed on Phase I to include 
rigid and flexure pivot tilt-pad (hydrostatic) bearings, as well as bearing shells 
with a simple elastic structural matrix (foil-bearings). The current THD model 
must be advanced to account for radial heat transfer through the bearing walls. 

There is also the need to develop a sound analysis and efficient computer 
code to calculate the transient fluid film bearing force response due to time- 
varying journal displacements and velocities. The analysis should be coupled to 
the dynamics of a rotor-bearing system for prediction of rotor lift-off, calculation 
of synchronous unbalance frequency response curves, and transient responses due 
to abnormal shock and maneuvering loads. 

The mechanics of two-phase flow phenomena and choked flow operation in 
unique cryogenic conditions should be addressed in the near future. A fundamental 
analysis is required with emphasis on the calculation of dynamic force coefficients 
and determination of stable regimes of operation. 




103 


references 

Of Hydrostatic JounSlfc^* Coeffidentf ” p 92 ’ “ E ^ perimentaI Determination 
International Conference in Vihrati • d r ? C * of IraecilE . Paper C432/145, 
September. ° ns ln Eotat ing Machinery, IMechE 1992-6, 

Arauz, G., and L. San Andres IQ*}*? “Pff * r 

S. .he Fee Response of a Short 

Journal Bearings with of H -'' brid Fluid Film 

Aided Bearing Design , ASME Press. ^ Effects ’ Adv ™ces in Computer 

S Fiow in Annular 

Fiow in Annuiar 

John wty & S So^«w 'Yo;k a ” d Lightfoot - E - N -. 19 «°. Tmnsporf Phenomena, 

VariaWe’prop^fc' ThemohyLu"!!: 198?a ’ “ A FuUy Cou P led 

Bearing,” ASME Joun„ AL OP 

™ 7 ?’ “Thermal Shaft Ef. 

Journal Bearing,” STLE Tribology Tmnlctioivl 30 “^.2^2-292” Cry ° g “ iC 

Pocket Aspect’ Ra«o, mdX SMph’cWfi’’ p"v " The E<feC ‘ S ° f a H y dr °static 
State FLow Patterns, Pressure .Id Shea rT a " d A “ act An « le °" Steady- 
Tribolooy, Vol. 115, pp. d Shear Ch “ueteristics,” ASME Jouanat. or 

Burton, R. A., and Carper H J 1 Qfi 7 u a i? 
with Applications in Turbulent Film LubrkaW ° f An T nular FI ° WS 

Technology, Vol. 89, pp. 381-391. ’ Journo ° f Lubrication 

a. Pump-Impeller- Shroud 

Stress and Reliabilify in ^eSgn Vo! ^ 1 ’!, Vibra,i °"' A ““«-s, 

Chaomleffel, J.P M and Nicholas, D. 1986 U F Y • . 

*’ 5 Experimental Investigation of 



104 


Hybrid Journal Bearings,” Tribology International, Vol. 19, No. 5, PP. 253-259. 

Constantinescu,V.N., 1973, “Basic Relationships in Turbulent Lubrication and 
Their Extension to Include Thermal Effects,” ASME Journal of Lubrication 
Iechnology, Vol. 95, pp. 147-154. 


Constantinescu,V.N., 
Force in Step Bearings,” 


and Galetuse,S., 1975, “Pressure Drop due to Inertia 
ASME Papers, 75-Lub-34. 


Di Pasquantonio,F., and Sala,P, 1984, “Influence of Thermal Field on the 
Resistance Law in Turbulent Bearing-Lubrication Theory,” ASME Journal of 
Iribology, Vol. 106, pp.368-376. 


Et ties , C . , and Cameron, A., 1968, “Considerations of Flow Across a Bearing 
Uroove, ASME Journal of Lubrication Technology, Vol. 90, pp. 312-319. 

Franchek,N., Childs, D., 1993, “Experimental Test Results for Four High-Speed 
High-Pressure, Orifice-Compensated Hybrid Bearings,” ASME Paper 93-Trib-24.’ 

Ghosh, M K Majumdar,B.C., and Rao,J.S., 1979, “Steady-State and Dynamic 
Behavim of^Mul^-Reces 5 Hybrid Oil Journal Bearings,” J. Mech. Engng. Set., 


Ghosh, M.IG, and Viswanath,N.S., 1987, “Recess Volume Fluid Compressibility 
uo^t 011 the i D y narmc Characteristics of Multirecess Hydrostatic Journal Bearings 
With Journal Rotation,” ASME Journal of Tribology, Vol. 109, pp. 417-426. 

Heshmat H., 1991, “Investigation of Foil Bearings for Use in High-Thrust 

Liquid Rocket Engines,” NASA Contractor Report CR-197099, MTI 90-TR58 
April. ’ 


Hill, D. E Baskharone, and L. San Andres, 1993, “Inertia Effects in a Hybrid 
Bearing with a 45 degree Entrance Region,” submitted to ASME Journal of 
Iribology, October. 

acSt??’ 0 ’’ 1973 ’ A Bulk ' Flow Theory for Turbulence in Lubricant Films,” 
ASME Journal of Lubrication Technology, Vol. 95, pp. 137-146. 

Ho,Y.S., and Chen,N.N., 1984, “Pressure Distribution in a Six-Pocket Hydro- 
static Journal Bearing,” Wear, Vol. 98, pp. 89-100. 

Holman, J.P., 1986, Heat Transfer (Sixth Edition), McGraw-Hill Book Company, 
New York. 

r o,^ 0nS ^V M< ’ 1987, “ A ReVieW ° f Therma] Effects in Hydrodynamic Bearings, 
1: blider/Thrust Bearings, II: Journal Bearings” STLE Tribology Transactions , 
Vol. 30, pp. 26-33. 

Knight, J. and A.J. Niewiarowski, 1990, “Effect of Two Film Rupture Models 


105 


Vol! n 2 Th p™M- A 188. ySiS ° f “ J ° UmiJ Beari " S ”’ ASME Jom,<AL °" Tr.BOLOGY, , 
Childs, D., San Andres,!., and Hale K 1993 “F T n • , v 

Z hZT Cha T eriStk : ° f a High - Sp - d Hybrid (Colina^r^dti,^ 

168. ynamiC ) Bearln 6) ASME Journal op Tribology, Vol. 115, pp.160- 

sSSSSsStr 21 “ “ S ■— 

College Station, TX 77843. ’ M y ’ Texas A&M Bniversity, 

r- a ~: *3* - - w 

Lubrication Technology, , Vol. 106, pp.237-245 ’ W 0r 

Va " “ “*“< < UK > Ca - 

Pmp^fe ’otnddf S MMOPS 86 f ” r Th e Dala Base * 2 ' 1986 ’ “Thermophysical 
Engineering, ^ C “*" Chemical 

Cooling" Effec!°j^ D a t d , 1988 ’ “Thermohydrodynamic Analysis of 

t ° bu PP ] y 0*1 in Circular Journal Bearing ” aqtup t 

Lubrication Technology, Vol. 105, pp.414-42™ *’ AL ° P 

Mori, A., Makino.T., and Mori IT 1 QQi U T?«+ 1 ^ 

Submerged Multi Part • j ’ 1991 ’, Entr y Flow and Pressure Jump in 
g Pad Bearings and Grooved Bearings,” ASME Paper 90-Trib-42. 

Mosher R, 1993, “Effect of Design Parameter Variations on Hybrid (Combi 

^L^Miter^h Hydr ,° dyn r C) BearingS fOT U “ High-Speed Turbo! 
77843. ’ eS1S ’ ay ’ Texas A&M University, College Station, TX 

with an' 4 ' ^ir^j’ ant * ^ Genge, 1993, “Hotfire Testing of a SSME HPOTP 

° f ^ ™ Workshop" 

A&M University, May. “ H ‘ gh P ' rformance Turbomachinery, Texas 

C —‘ Heat,” Meehan- 


106 


Pinkus,0., 1990, Thermal Aspects of Fluid Film Tribology , ASME Press, New 
York. 

Reddecliff,J.M., and Vohr,J.H., 1969, “Hydrostatic Bearings for Cryogenic 
Rocket Engine Turbopumps,” ASME J ournal of Lubrication Technology, 
Vol. 91, pp. 557-575. 

Roach, P.J., 1976, Computational Fluid Dynamics , Hermosa Publishers, Albu- 
querque, N.M. 

Rohde, S.M., and Ezzat,H.A., 1976, “On the Dynamic Behavior of Hybrid 
Journal Bearings,” ASME Journal of Lubrication Technology, Vol. 98, pp.557- 
575. 

Safar, Z., and A. Z. Szeri, “Thermohydrodynamic Lubrication in Laminar and 
Turbulent Regimes,” ASME Journal of Lubrication Technology, Vol. 96, pp.48— 
57. 

San Andres, L., 1990a, “Turbulent Hybrid Bearings with Fluid Inertia Effects,” 
ASME Journal of Tribology, Vol. 112, pp. 699-707. 

San Andres, L., 1990b, “Approximate Analysis of Turbulent Hybrid Bearings: 
Static and Dynamic Performance for Centered Operation,” ASME Journal of 
Tribology, Vol. 112, pp. 692-698. 

San Andres, L., 1990c, “Improved Analysis of High Speed, Turbulent Hybrid 
Bearings” 4th NASA Conference on Advanced Earth-to-Orbit Propulsion Tech- 
nology, NASA CP 3092, Vol. 2, pp. 414-431, May 

San Andres, Luis A., 1991a, “Effect of Fluid Compressibility on the Dynamic 
Response of Hydrostatic Journal Bearings,” Wear, Vol. 146, pp. 269-283. 

San Andres, Luis A., 1991b, “Effect of Eccentricity on the Force Response of a 
Hybrid Bearing STLE Tribology Transactions, Vol. 34, 4, pp. 537-544. 

San Andres, Luis A., 1991c, “Analysis of Variable Fluid Properties, Turbulent 
Annular Seals,” ASME Journal of Tribology, Vol. 113, pp.694-702. 

San Andres, Luis A., 1992a, “Analysis of Turbulent Hydrostatic Bearings with 
a Barotropic Cryogenic Fluid,” ASME Journal of Tribology, Vol. 114, pp.755— 
765. 

San Andres, L., 1992b, “Analysis of Hydrostatic Bearings with End Seals,” 
ASME Journal of Tribology, Vol. 114, pp. 802-811. 

San Andres, L., 1993a, “The Effect of Journal Misalignment on the Operation of 
Turbulent Flow Hydrostatic Bearings,” ASME Journal of Tribology, Vol. 115, 
pp. 355-363. 

San Andres, L., 1993b, “Effect of Shaft Misalignment on the Dynamic Response 



107 


of Annular Pressure Seals,” STLE Tribology Transactions , Vol. 36, pp. 173-182. 

San Andres, L., 1993c, “Dynamic Force and Moment Coefficients for Short 
Length Annular Seals,” ASME Journal of Tribology, Vol. 115, pp. 134-140. 

San Andres, L., 1993d, “hydrosealt User’s Manual and Tutorial,” Texas A&M 
University, Mechanical Engineering Department, December. 

San Andres, L., 1993e, “hydrosealt Examples,” Texas A&M University, Mechan- 
ical Engineering Department, December. 

San^ Andres, L., 1993f, “Turbulent Flow Foil Bearings for Cryogenic Applica- 
tions, submitted for review to ASME Journal of Tribology, December. 

San Andres, L., and Velthuis,J.F.M., 1992, “Laminar Flow in a Recess of a 
Hydrostatic Bearing,” Tribology Transactions , Vol. 35, pp. 736-744. 

Scharrer,J.K., and Hibbs,R.I., 1990, “Flow Coefficients for the Orifice of a 
Hydrostatic Bearing,” STLE Tribology Transactions , Vol. 33, pp. 543-550. 

Scharrer,J.K., R. Hecht, and Hibbs,R.I., 1991a, “The Effects of Wear on 
the Rotordynamic Coefficients of a Hydrostatic Bearing,” ASME Journal of 
Tribology, , Vol. 113, pp. 210-21. 

Scharrer,J.K., J. Tellier, and Hibbs, R.I., 1991b, “A Study of the Transient 
Performance of Hydrostatic Journal Bearings, I: Test Apparatus and Facility, II: 
Experimental Results, ” STLE Tribology Transactions , Preprint 91-TC-3B-2 

Scharrer,J.K., and T.W. Henderson, 1992a, “Hydrostatic Bearing Selection for 
the STME Hydrogen Turbopump,” AIAA Paper 92-3283, 28th AIA A/S AE/ ASME 
/ASEE Joint Propulsion Conference, July 6-8, Tennese. 

Scharrer,J.K., Hibbs, R. I., and San Andres, L., 1992b, “The Axisymmetrically 
Stepped, Orifice Compensated Hydrostatic Bearings,” AIAA Paper 92-3405, 28th 
AIA A/S AE/ ASME/ ASEE Joint Propulsion Conference, July 6-8, Tennese. 

ce?™ a £™’ K '’ Hibbs, R.I., and Nolan, S., 1992c, “Extending the Life of the 
SSME HPOTP through the Use of Annular Hydrostatic Bearings,” AIAA Paper 

92-3401, 28th AIAA/SAE/ASME /ASEE Joint Propulsion Conference, July 6-8, 
Tennese. 

Scharrer,J.K., J. Tellier, and R. Hibbs, 1992d, “Start Transient of an An- 
nular Hydrostatic Bearing in Liquid Oxygen,” AIAA Paper 92-3404, 28th 
AIAA/SAE/ASME /ASEE Joint Propulsion Conference, July 6-8, Tennese. 

Shoup, T., Private Communication, 1993, Sverdrup-NASA Marshall Space 
Center, October 15. 

Simon, J., and J. Frene, 1989, Static and Dynamic Characteristics of Turbulent 



108 


iZ'pr ECC ! n,r » C A S ™ , E<feCt ° f Gonver 8 en I 'Tapered Geometry and Variable 
r luid Properties” ASME Journal of Tribology, Vol. Ill, pp.378-385. 

Smalley, A.J., Vohr,J.H., CastelU.V. and Wachtmann,C., 1974, “An Analytical 
and Experimental Investigation of Turbulent Flow in Bearing Films Including 

Vol 96 po 15 l‘t 5 7 ert ‘ a ASME JoUnNAl or Technology., 

Suganami.T. and Szeri.A-Z., 1979, “A Thermohydrodynamic Analysis of Jour- 
nal Bearings, ASME Journal of Lubrication Technology, Vol. 101, pp.21-27. 

Tonnesen, J. and P K. Hansen, 1981, “Some Experiments on the Steady State 

° f “ 5 y ‘ ndnCal Fluid -Film Bearing Considering Thermal Effects,” 
oi Ul Journal of Lubrication Technology, Vol. 103, pp. 107-114. 

MTr J D r 0r o aal ,’ J ' R ’ and Raithb y. G D -, 1984, “Enhancements of the SIMPLE 
Method for Pred.ct.ng Incompressible Fluid Flows,” Numerical Heat Transfer, 
Vol. 7, pp. 147-163. 

Van DoormaaU.P., and Raithby.G.D., 1985, “An Evaluation of the Segregated 
Approach for Predicting Incompressible Fluid Flow,” ASME Papers, 85-HT-9. 

19^ Pra6enau,G-L -’ 1982 ’ “Damping seals for Turbomachinery,” NASA TP 

Qnoo n Af[ a ?? naU ’ G ' L ‘’ 199 °’ “ Dam P in g bearings for Turbomachinery,” NASA TP 
3092, Vol. II, pp . 155-162. 

In c W NlwYfrk 1977 ’ FUTldamental Princ 'P l ™ of Heat Transfer , Pergamon Press 

Yang Z., San Andres, L., and Childs, D., 1993a, “Thermal Effects in Cryogenic 

• 19l “ d Annular 8eals ’ Part I: Theory and Approximate Solution; Part II: Numer- 
ical Solution and Results” ASME Journal of Tribology, Vol. 115, pp.267-284. 

Yang Z San Andres, L., and Childs, D., 1993b, “Importance of Heat Transfer 
from Fluid Film to Stator in Turbulent Annular Seals,” Wear, Vol. 160, pp.269- 

Ying,Z., San Andres, L., and Childs, D., 1993c, “Thermohydrodynamic Analysis 
of Process Liquid Hydrostatic Bearings in Turbulent Regime, I: Theory, II: 
umencal solution and Results”, accepted for publication at ASME Journal of 
Applied Mechanics, December. 

Yang Z San Andres, L, and Childs, D., 1993d, “Process Liquid Turbulent Flow 
HydrosUtm Bearings: Analysis and Tests for Rotordynamic Coefficients”, ac- 
cepted for presentation at 1994 IFToMM Rotordynamics Conference, December. 

Zeidan, F., 1992, “Developments in Fluid Film Bearing Technology,” Turboma- 


chinery International, September/October. 




110 


APPENDIX A 
NOMENCLATURE 


A 

itDL, journal or bearing surface area [ra 2 ] 

A 0 

Cdirdl/4, equivalent orifice area [m 2 ] 

A r 

bl, recess area [m 2 ] 

b 

recess circumferential length [m] 

c^c* ,c 

radial clearance, characteristic clearance (={c(y)}mm) [m], c/c. 

c d 

empirical orifice discharge coefficient 

c p ,c p 

specific heat [J/kg -K},C p /C pm nondimensional specific heat 

^ 'xx 

damping coefficients 

D 

journal diameter [m] 

d 0 

orifice diameter [m] 

e *><Y 

displacements of the journal in X and Y direction [m] 

E c 

U? /T*C pm , Eckert number 

f:J B 

a m [l + {c m r J B /H + turbulent friction factors at 

journal and bushing surfaces based on Moody’s equation, 
a m =0.001375; 6 m =5xl0 5 ; c m =10 4 ; e m = l/2.65 

frj 

Sr. h <“> 

F x ,F y ,M x ,M y 

film forces(moments) along(around) {X,Y} ajces [N,N-m] 

H , h 

film thickness [m], H /c* 

Hr 

recess depth [mj 

ho 5 h i j 

zeroth- and first-order film thicknesses 

hg,hj 

heat transfer coefficients to bearing, journal surfaces [W /in 2 • K\ 

h B J 1 ] 

dimensionless bulk heat transfer coefficients 

h x ihy ih$ x 

cos 0, sin 0, sin 6, cos 6, perturbed film thickness components 



Ill 


K,K 

^ XX T'-iKSydy 

k k 

a I jA/y 

k iA 

k r 

Li LliLr 

l 

m 

M 

M r 

Kxi-iM^ 

n 

Pip 

PxiPy 

Pe 

Pl,Pr 

PsiPaiPr 


thermal conductivity of liquid [W/m -K), K/K, 
stiffness coefficients 


(K + k B )/2, dimensionless shear parameters in X, Y directions 

fjR ^' tUrbulent sbear parameters at journal, bearing surfaces 
Pe°r- 661 / 7 . 753 , turbulent shear flow parameter 
bearing axial length [mj, L = L l + L R 
recess axial length [m] 

mass flow rate over differential segments [kg/s] 


mass flow rate of a hydrostatic journal bearing [kg/s] 

recess flow Mach number due to rotation 
inertia or added mass coefficients [kg] 

normal vector to recess boundary 

fluid pressure (p _ P o)/Ai >, dimensionless pressure 


imens.onless dynamic pressures for perturbation 
te‘vf>r = P-U.clCp./(K.R), Peclet number 
discharge pressures on left, right sides of bearing [N/m 2 ] 
external supply, ambient and recess pressures [N/m 2 ] 
C p fi/K, Prandtl number 


mass flow rate supplied to a recess/groove [kg/s] 

Qin/{p*C.RU*) 

recess to land mass flow rate [kg/s] 

Qr/(p*c m RU *) = J ph( u • n)dr r 

heat flux to the bounding surfaces [W/m 2 ] 

dimensionless heat flux to the bounding surfaces 
journal radius [m] 

P.Wc./p., nominal circumferential flow Reynolds number 
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Re H 

PeRftHe/ve, Reynolds number at recess edges 

Re p 

p.U,c./n., reference pressure flow Reynolds number 

Re* p 

Re p (c*/R ), modified reference pressure flow Reynolds number 

Re r 

PrRfl(H r + H)j Hr, Reynolds No. at recess volume due to rotation 

Re, 

fi,, nominal squeeze film Reynolds number 

Rj 

pHy/(U - QR ) 2 + V 2 /fi Reynolds number to journal surface 

Rb 

PH'/U 2 + V 2 /p, Reynolds number relative to bearing surface 

r j’ r B 

mean roughness depth at journal and bearing surfaces [m] 

S t 

ht/pCpVt, Stanton number 

T,T 

two-dimensional bulk temperature in the fluid film [A'], T/T » 

T. 

T si characteristic temperature [A"] 

T 

T{x,y,z,t), three-dimensional temperature [A"] 

t 

time [sec] 

Tb,Tj 

temperatures at the bearing and the journal surfaces [K] 

T, 

inlet supply temperature at the orifice [K] 

Tu ,T d ,T 3lde 

temperatures at the up-, down-stream and side recess edges [K] 

T or 

fluid film resistant torque [N-m] 

rpv 

-*■ or 

T xzA r R, friction torque over a recess area [N-m] 

u,v 

mean flow circumferential and axial velocities [m/s] 

u,v 

(U>V)/£/*, dimensionless mean flow velocities 

U,Vr 

Ui + Vj, RCli + 0/ 

u,v,w 

three-dimensional velocities in x,y,z directions [m/s] 

u. 

^R c */ (p* R), characteristic velocity [m/s] 

V 

{U,V,W} 

V. 

(Rr + H'jAr + V„ total recess volume [m 3 ] 

Vr 

V r /(c*A 2 ), dimensionless recess volume 
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V, 

V t ,v t 

X,Y,Z 


x,y,z 

ZoX,Z 


x,y 

fi P 

fi P 

fit 

fit 


a 

a y 

A 

<7 

A 

Ai 

^11 

A2 

n 

u? 

P)P*,P 


Pip* 

<t>r 

4>x,y 


volume of orifice supply line [m 3 ] 

y/U 2 + V 2 [m/s], \/u 2 4- v 2 , fluid speed 

inertial coordinates defining journal position in bearing 

(0,7rD),(~£^,£ fi ) j (0,H(x,y,t)), coordinates defining flow regions 

distance to center of journal angular motions, Zq/R , Z/R 

( X /R = Q),ylR, dimensionless x, y coordinates 

p(dp/&P) T i liquid compressibility factor [m 2 fN] 

dimensionless liquid compressibility factor [m 2 /N] 
fiT) P i volumetric expansion coefficient [ 1 /A'j 
fitT», dimensionless volumetric expansion coefficient 
( ^ | j/=o )/(QR), circumferential velocity entrance swirl factor 
2(1 ~b £)Re p , axial entrance loss coefficient 
£IR/U., dimensionless journal velocity or speed parameter 
uiR/U*, dimensionless frequency parameter at film land 
thermal mixing coefficient for groove heat carry-over 
crV r , dimensionless frequency parameter at recess volume 
crl/ R, modified frequency parameter at recess volume 
^ 1 fip r ) combination of compressibility & frequency parameters 
rotational speed of journal [rad/sec] 
excitation or whirling frequency [rad/sec] 
fluid density, characteristic density [kg/m 3 ], p/p , 
fluid viscosity, characteristic viscosity [Ns/m 2 ], p/p* 

recess volume shear coefficient 

^ e p(l + Cc,y)(l + £z,y)[l ~ (vPe/ pr) 2 ], recess entrance loss factors 
1.95 /Re^\ 0 

empirical entrance loss coefficients in X, Y directions 
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€xu,d 

6* 

e x 

KA 

Ae x’^ 

^xAS, 

r 

r c k 

r xz,T yz 

7 

7 

K P 

KlXr 

Scripts: 

0 

3 

e 

r 

J 

B 

* 


a ‘ “ P ' and do "'"- s ‘^™ of recess/land entrance, respectively 
A^.s/2/l ciVp^Fl, dimensionless orifice parameter 

dimensionless jonrnal eccentricities in X, Y directions 
journal axis angular rotations around X, Y axes 

dimensionless dynamic (perturbed) eccentricities 
dynamic (perturbed) journal angular rotations 
wt, dimensionless time coordinate 

riuk B ~( u ~ A)k J j/4h, dimensionless shear stress on journal 
wall shear stresses in X and Y directions 

H/(H r + H), ratio of land film thickness to recess depth 
first order turbulent shear coefficients 

pad leading edge pressure recovery factor 
discharge coefficients for end seal restrictions 

refers to zeroth-order solution 

refers to first-order perturbations), - X. V direction) 

refers to entrance or recess edge conditions 

refers to recess conditions 

refers to journal 

refers to bushing 

refers to characteristic values or supply conditions 

refers to three-dimensional variables 

refers to dimensionless variables or parameters 
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appendix b 

heat transfer to bushing and 

For the lumped fluid film model 


JOURNAL 


assumption, the heat transfer from fluid film 
to the bounding surfaces has been expressed 


as 


Q, = hB(T-T B ) + hj(T-Tj) (B1) 

^ “ d ^ ale tHe coefficients to the bushing and journal 

respectively. They can he found through the Reynolds-Colburn analogy between 

flu.d fr.ct.on and hea, transfer (Hohnan, 1986). The average hea, transfer over 
the entire laminar/turbulent boundary layer is 


where 


StP 2 J 3 = f/2 


Q *lt 

* ~ pC p V t ( Stanton number ), 

_ CpfX 

^ T ~ ft ( Prandtl number ), and 

'-M 1 + («-£ + £)"] 

is the Fanning friction factor Based on Moody’s Equation (Massey, 1983). 
So, the heat transfer coefficient can be written as 

h < = \ec,v,nvy' 

= \pC P v B f B if>y a 


(B. 2) 

(B. 3) 
{BA) 
(B.S) 


(B.6) 
(B. 7) 


* = 2 Wh/A* (B.8) 

Actually, there are plenty of formulae for turbulent heat transfer coefficients 

(Holman, 1986). The above expressions are used because of the simplicity and 
the ab.l.ty to include the surface roughness effects. 
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APPENDIX C 

ZEROTH- AND FIRST-ORDER EQUATIONS 
C.l Zeroth-Order Equations on Film Lands 

• Continuity Equation 


9 Q 

-QzlflohoUo) + -Qz(poh Q v 0 ) = 0 ; 


Circumfere ntial- Momentum Equation 

d Po _ , A, „ 9 


-h 


°~dt = T^ kM ~ kr0 2^ + Re p\-lk(P* h * u o) + nzipohouovo )} ; (C.2) 


dy 


Axial-Momentum Equation 


$j) - 

h °~dy = D ky0V ^ + Re *pifc(pohou 0 v 0 ) + ^z(p 0 h 0 vl ()] ; 


dy u 0 

• Energy- Transport Equation 

Re P f n r d / -9 i 

~E7\ C Poig^(po h 0 u 0 To) + Qz(pohov 0 T 0 )} + ( h Bo + A Jo )f 0 J 

3 rp L / dPo dpn . A dp n R^n - - - 

- PtoT 0 ho{u 0 — + v 0 ^r-) + ho--? 1 + — ^(h T +h T ) 
ox dy 2 dx E c ' 80 s» T V;oi 

fl 0 i 2 

+ — uoA) + A: r o( — A 2 — uqA)] 

C.2 Zeroth-Order Equations at HJB Recesses 

• Mass Conservation at a Recess 


(C.l) 


(C. 3) 


(CA) 


fi*V~Pro(l — Pro ) — 9r0 = / /?0 A 0 (u 0 • n) dF 

Jr r 

• Global Energy BalanrP Equation at a Recess 
Cpo (^^ di^Tda + 

Inside 0 ^aideO^j 

=Cpo (XI ™uof uQ + 9tn0 ) + (f; r0 A)E c /Re; 


(C. 5) 


(C.6) 
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The recess temperature can be found from the above energy balance equation 
as follows. 

In a hydrostatic journal bearing, flow supply enters the recess through an 
orifice restrictor. Temperatures at the downstream and side edges of the recess are 
approximately equal to the recess temperature due to the high pressure gradient: 


Tdo — T s ifi e o — T r0 = Constant 


(c. 7) 


Substituting Eq.(C.7) into Eq.(C.6), the recess temperature is then given by: 

1 


TrO — 


0 n ~ [?m0 + y^m u0 T u0 + 

2<7^eo + <ld0 L ^ ^ C p o Re; J 


(C.8) 


where 


T u q 


Tr 0 , 

Upstream Film Values on Film Lands, 


if (U 0 • ii) > 0; 
otherwise. 


(C- 9) 


The above equations show that if fluid flows out of all the four recess edges 
(as in a pure hydrostatic beairng), the edge temperatures are equal to the recess 
temperature. If the hot lubricant left the upstream film lands enters the recess 
at the trailing edge due to journal rotation, then an averaged uniform recess 
temperature can be obtained from the mixing of upstream and inlet temperatures. 

# Global Energy Balanc e Equation at a fironvp 


where 


If 


zlw 7 ^ ld oT'dO ^ ] ri^uoTuO + qinoTt'nO 

(C.10) 

?rn° = ^2 ™<io - ^ m u0 = q d0 - q u0 . 

(C.ll) 

Quo 1 ^ m u0 > o, 

(c. 12) 


then, 
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2uo Upstream Film Values on Upstream Pad (C.13) 


and 


( a ) U 9rfo > o, then 


v qdo * 


9d 0 


inO 


where 


T in0 = 1 . 0 . 


(C. 14) 


(C.15) 


(b) If 

<7<*o < <7uo> then 


f XjT u q^ .£ ^>0* 

1 Upstream Film Values on Downstream Pad, if ^ < 0.' (C 1 - 16 ) 


* Recess /Film Entrance Pressure Rise/Drop 
Peo = I Pr0 ~ ^rotrA ' - y] > (u 0 • n) > 0, in x direction 


PrO, 


in y direction; 




•n ' _ & * PeO , - 2 

Peo-Peo g ( u o- n ) , ( U 0 • n) > 0 


(C.17) 


(C.18) 


C.3 First-Order Equations on Film Lands 

The fist-order equations are obtained from the perturbations in the film 
thickness due to journal displacements (Ae x ,Ae r ) and rotations (A6 X ,A6 X ). The 
perturbed (dynamic) film thickness in the following first-order equations can be 


expressed as 

h 


. _ / KA’ 

‘ ~ \ (z - z„ 


f° r journal displacements about (X,Y) axes; 
){ h s x , hsy }, for journal axis rotations around (X,Y). 


where 


(C.19) 


h x = h 6y = COS 6, hy = -h s = sin0 


{C. 20 ) 
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• Continuity Equation 

d _ q 

gzipouohj + Poujho + pjuoho) + Qzipovthi +p 9 Vjh 9 ^p jV ^) 

+ ia(p 0 hj + pjho) — 0 

• Circu mferential-Momentum Equation 
-h ^2l _ 

°~dF - ' ! ' ,,L + + 1,,Vj + 7,A Aj + JvpPj + luppj 


(C. 21 ) 


+ + A(M.»ou>) + pohoVj <C ' 22) 

• Axial-Momentum Equation 


, C>Pj 

= (7uv + i Mofle,)u i + 7.t,tt; + 


pPj ' IfvfiPj 


Q Q 

+ Re *p{fa(f> oh <> u <> v j) + Qg(pohoV 0 Vj) + Pohou^ + p 0 hoVj^-] 
Energy - Transport Equation 
Re* Q q 

C p °~E7 ^P° h o u oTj) + gzipohovofj)} + ( 7tt 4- ip 0 h 0 C Po Re./E e )Tj 

+7tuUj + 7 tvvj + ( 7 t P - i<T0 to foho) Pj + jtkhj - ^ 0 - — 

2 dx 


(C. 23) 


(C.24) 


- PtoT 0 ho(uo^r + v 0 ^-) + T + h T ) 

here, the subscript j (— X ,Y,6 X ,6 y ) denotes the direction of perturbation for the 
first-order variables. The 7 coefficients arise from the perturbation of the turbulent 
shear parameters k x and k y and are given in this Appendix. The perturbed first- 
order liquid properties are evaluated from relationships like 




d<j>] 

Idp J 

o Pj + 

IdrJ 


Tj , 4> = p,p,...,etc. 


(C. 25) 


C.4 First-Order Equations at HJB Recesses 
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• Mass Conservation at a Recess 


prj | 2 [l ^pro(l — PrO )] + 1^2/^rO 


r fi 2 p o 

z Qrj “h i^rO^ll frj + (l — Pt-o) - ^ l^lPrO\0troTrj 

L lq r o J 


(C. 26) 


where 


■j — I hj 

J r. 


<j r j — / (po^j u o + pohgiij -(- pj H q Uo ) ” n , 

J r. 

An = crl/R , 

A 2 — Ai/?p r o — frV r( t3p r o 


(C.27) 

(C.28) 

(C.29) 

(C.30) 


are the perturbed recess mass flow rate, a modified frequency parameter, and a 
combination of the frequency and compressibility parameters at recess volume, 
respectively. 

• Global Energy Balance Equation at a Recess 

^ d 0 Rdj ”1" TTldj'RdO 7Tl u Q r L' u j ^uO d" ^^TlsideO^sidej d~ *^771 s ide sideO^j 

_ i Q /'T - i — 'v \ i f'^' or i Cprj ^ T or0 E c A 

— q r j UrV r (T r op r j + p r oTrj ) + ( ^ ^ ^ 77 ^ 

Vi or -0 L 'prO / (-'prO ne p 


Assuming 

^uj = Tdj = = T rJ , 

then, the perturbed recess temperature can be calculated as 

'^'rj . q I 9rj ‘2q a idejT r Q + ^^^T7l u jT u 0 ^^^TTldjTdO 


((7.31) 

(C.32) 


9r0 + i(7V r p r o 


rpr 

( or j 

Cprj N 

\ 2 or 0 

EcAi 

\T r 

W or 0 

o 

^ prO 

;? e ;J 


(C.33) 
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where 




7 s^m \ Pj 
ho ho 


/ 1 \ pj / 2 Cca \ hj j 

' ho ho ' ho h 0 \ r 


is the perturbed dimensionless torque. 


• Global Energy Balance Equation at a Groove 


(C. 34) 


Assuming 

( T in )j = 0 

then 


y: ( rhdoTdj + rh d jT d 0 - m u0 T u j 


riiujTu o 


) = £(" 


riiuj ) Tino ((7.35) 


if 


9u0 


= E 


m u0 > 0, 


then 

Tuj = Upstream Film Values on Upstream Pad, 

and 

(a) If q d o = Y, ^do > q u o , then 


(C.36) 


Tdj — — ^ A (rh u0 T u j + rh u jT u0 — rh d jT d 0 ) + (q d j - Xq u j)T in0 (C.37) 


(b) If q d o < q u0 , then 


2^. _ J ATuj, if > 0; 

\ Upstream Film Values on Downstream Pad, if q d 0 < 0. 


((7.38) 


• Recess/Film Entrance Pressure Rise/Drop 


P 


ej 


{ 


p rj - Pr<t>rV{j;)(Uj • fi), 
Prj j 


Uo * n > 0, in x direction 
in y direction; 


(6739) 
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Pej Pej <f>i{peo(u 0 ■ n)(fij • n) + p e j(u 0 • n) 2 /2] ; u 0 • n > 0 (C.40) 


C.5 Perturbed Shear Coefficients 


The perturbed shear coefficients for the first-order solution are 

given for a 

general fluid by: 



f ^ 21 0 n 

7,i “ = f + ^ rl ( u ° _ A ) + 

(C.41) 


7v " ~ + ^ rl + 

(<7.42) 

7ut> — — [/rl(«0 ~ A) -f /,i U 0 ]t’0^0 

((7.43) 

luh — 

“«(“** o 7- C cr + C cs ) - A(C cr - * J0 /2)] 0- 

((7.44) 


^0 

7u/i = 7y( — h y o + C cr -)- <7 c ,)/io 
n o 

(C.45) 

7 up 

= kr^ + ->- 

^"0 r0 /? 0 C7x 

(<7. 46) 


Ao r , X1 ho dfl, 

7** 2 hn n M 7 .+ 73 )] _ 

Po /?o c/y 

((7.47) 


r , 

TV ~ 2h 0 ^ U °^ r - AX7,] 

((7.48) 


Twm — 2^ 0 f U °(7r- + 73 )] 

((7.49) 

7 ‘“ = + + - A >,?j} 

BO K j o ^ 



If tv 


Re; 

' E c 


^ w j u 

- 0toToh o ■— — El[ k x0 (2u 0 + y) - AA^ 0 ] 
^o[(«co + 2u c i)(u 0 — A )/ rl + u C QU(j f 3 i ] 

df ° -• '/•!* . frl ^ ^ ^ 


< 9 y 




y 0 

~ -j^( 2k yoVo) -/ioVo[(«co(/n + fsl) + 2u cl f rl ] 

?P . „ d fi , dK dc p 

1,p - ^ a? + ^ ™ af + + t c , 


(<7.50) 


(C.51) 


(<7. 52) 
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Re' 

T " = £T (/! *" + ^-A *.(«.&+*&) 

d * 9y> 

+ 7<»jj7 + 7l + 5/j ea ac 

dt 7 <« g, + 10 -$■ + 7tt ~ + lc f^p 


(<7.53) 


dt lc v fa 

7 ' 4 + "*f > + £( §«.. + 

- - <9n ,9 so v° ■* 2 <9# 

■PoT 0 (u 0 ^.+ Vo %_h' (C. 54) 

0* 0 <9v ' Ty l^co(O c > S' ' ■ ~ - - ' 


+ <? e.)+2u cl C er ] 

*. ___ r j- ^ 

2 Uc0 ^ 7r + 7i ) + u ci7r] ~ -z~~ f-Q J0 + + 7rbm n 1 

MoR c l3 Va0+ k Bo Q *° + ~QjoJ (C. 55) 


di °dy } hi 

bmrl ' r2 


7t P 


~~ p -\r u / < 9 r 0 #f n 

L ?0 o(&0 ~ + l?n 1 _L /.Q , 7 ab m - -y ft 

f*o + — — <0 + w - i 

^0 fl0 Ar. ^o)/poj 


/aq f ^ fl0 

Mo <M **» + f + 7.)1 + ««(*,„ + 

TW = -A,1'.( 11 ,$| + 

(9a: Qy } 

7tk = 

3 £ c 


7»^ m _ 


/O 


iZe* 

~E, 


where 


= * l *" Cr + 6 ” hv ; c - = + 4mj7 „ 


C r = c m -III. . r , r s 

m r h ' C * ~ c m — 

Cmh ° C,h 0 ’ 

R jo = R e p(po/ p,o)h 0 f (^ 0 _ A ) 2 + v *ji/a 

= Re p(po/fio)h 0 [u 2 0 + ^j 1 / 2 , 

" a 7n * e m 


(C.56) 

(C.57) 
(C. 58) 
(O. 59) 

(C.60) 

(0.61) 

( 0.6 2 ) 


7 r = 
7s = 


[fro/a m - l](l/e m -l) ’ 

~ Q 7 n • e m 

[fso/dm — lj(l/ e m~l) * 


(<7.63) 
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f , _ (*WM 2 V, 

2f> [fr 0 + /?r7rj ] 



f 1 _ (^ e pPo//io) 2 /io r r . 

‘'* 1 2i^ [fsO+Psls], 

((7.64) 


O 

II 

•> 

o 

g 

~C> 

II 

((7.65) 


~ “» + ®J + “.f ; “«i = - u»A, 

((7.66) 

- 

= Q,. = Q„„ + Q J0 

(<7.67) 



125 


APPENDIX D 

ALGEBRAIC EQUATIONS FOR THIN FILM FLUID FLOWS 


For convenience, the subscript “0” is omitted for all the zeroth-order variables. 
D.l Dimensionless Zeroth-Order Discretization Equations 
a) Circumferential u-momentum equation : 


^P u p ~ ^ U E + + ^ s u s + ^ u n + ^ (D.l) 


where 


4 = !>::»+$? 

'£<i, = A U e + A*w + A'‘s + A’‘ n 
A% = fle;Max(-F e “,0); A" w = SeJMaxff^.O) j 
A' n = i?€pMax( - F" , 0) ; A‘ s = fl e ;Max(F",0) 

Fg — ^phu^j^Ay , Fy, — (phu)mAy ^ 

F n = (p/iv) n Ax ; Fj = (p/iv) s Ax 


= 


r tt k 1 

ih kx \ 


\Ax u Ay u 


5“ = 


tK- 

U 2 J p 


‘ Ax u Ay u - A£(p e - p w )Ay u 


b) Axial v-momentum equation : 


(D.2) 

(0.3) 

(0.4) 

(0.5) 

( 0 . 6 ) 

(0.7) 


4.» F 


+ Aw v w + ^s v s + 


(D. 8) 


where 


4 ;, = £ a; 4 + sj 


(0.9) 
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A v e = i2e*Max(— F" , 0) ; A v w = fl^Max^O) ; 
A v n = iZe*Max( — F",0) ; .4£ = i2e*Max(F, v , 0) 


Sl = 


t i k 
L h y 


v p Ax v Ay v 


s c = ~ h p(p n ~P,)Ax v 

c) Correction-p equationf from the continuity equation 

ApP P = A p E p e + A p w p w + A p s p' s + A p N p' N + S p 

A ? p = ^2 A^ b 

A p E = (ph) p Ay p D p ; A* = (ph) p Ax p D; ; 


where 


D“ 


^Ay p 




F>! = 


/tpAxp 


jt)_r ju ’ ' 
2^ ^nf> 


Note that the denominators of D? and D * will not be zero due to 
relaxation and source terms. 


s; = -(j? -fz+k- k) 


d) Energy equation : 

ApTp = A e T e + A* W T W + A^Tg + .4^7^ + S* 

A* p = ^2 + Spi + Max(5p 2 , 0) 

A e = C p (i2e;/F c )Max(-F e t ,0); A f w = C p ( J R e ;/F c )Max(F^ 0) ; 
a n = Cp(Re*/ E c )Max.(—F^, 0) ; A* s = C p (i2e;/F c )Max(F;,0) 
5pi = Re*/E c (h g + hJAipAj/p 

*^P2 — dp/ip[zip(p e — p w }Ay p + Vp{p n — p,)Ai p ] 


(O.10) 

(O.ll) 

(D.12) 

(0.13) 

(0.14) 

(0.15) 

(0.16) 

under- 

(0.17) 

( 0 . 18 ) 

(0.19) 

( 0 . 20 ) 

( 0 . 21 ) 


(D.22) 
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St = S l cl + 5* 2 + 5*3 + Max(— 5p, 0)Tp O) (£>.23) 

Sci = {^[kx(u 2 + y2 + 2^^^ + ^r( — uA)] j^AxpAy p (D.24) 

S c 2 ~2^p{Pe Pu;)Aj/p (Z). 25) 

5*3 = Re;/E c (h B f B + hjT^AxpAyp (D. 26) 

D-2 Dimensionless First-Order Discretization Equations 

a) u-momentum equation : 

A Pj u jP = A E iij E + AwUj W -l- A5 «. s + A^u^ + S*j (D. 27) 

where 

A Pj = lL A nb + S Pj {D.28) 

A\ = Z?e*Max( — , 0) ; = Z7 e ;Max(F“, 0) ; 

(D. 29) 

^)v = i2e*Max(-F“,0) ; A£ = A^Max^, 0) 

Sp; = (7uu + iphRe a )Ax u Ay u + Max(5£ 2j ,0) (D.30) 

^P 2 j = Re* p ph(u t - «u,)Ay u (D.31) 

Scj = - (luvVj + -fuh.hj -f TupPj + 'j U) ifi'j)Ax u Ay u 

— Re*pk(u n — u a )Ax u ir,- + Max( — 5p 2 j, Ojtij (D.32) 

— hp(pj E — p. p )Ay u 

b) v- momentum equation : 

A Pj V jP = A E v jE + + ^S u jS + + $cj (D- 33) 

*Pj = £ AJU + 5pj (£>.34) 


where 
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A E = «e;Max(-F e u ,0); A v w = fl e ;Max(F”, 0) ; 
A n = Re p Max( — F", 0) ; A v s = Fe*Max(F i t ’, 0) 

(£>. 35) 


Spj — (iw + iphRe a )Ax v Ay v + Max(5£ 2j , 0) 

(£>.36) 


Sp 2 j = Re* p ph(v n - v,)Ax v 

(£>.37) 


Scj = - (jvuUj + 7 vkhj -f 7 vpPj + ' Y Vfl p.j)Ax v Ay v 



— Re*ph(v e — v w )A y v Uj + Max( — 5p 2j , 0)uj 

- h v P (p jp - p jS ) Ax v 

(£>.38) 

Correction-p equation( from the continuity equation ): 



A pPjP A ePjE + A wPjW + A sPjs + A nPj n + S cj 

(£>.39) 

where 

A P = XI A nb 

(£>.40) 


A P e = (ph) p e Ay p D u p ; A p s = (ph) p Ax p D v p ; 

(£>.41) 


P . _ K*«r p , 

' A i-Y.Ki' p A l-Y.KP 

(£>.42) 


S cj = - ( F ej - F tj + F nj - F sj ) - (puhj + pjuh) e w Ay 
— ( pvhj + Pjvh)*Ax — i cr(pjh + phj)AxAy 

(£>.43) 


F ej = (phuj )e A y] F P J = ( phuj ) w Ay ; 
F nj = ( phvj) n Ax ; F P j = (phvj),Ax 

(£>.44) 

Energy equation: 



A Pj T jP - A E T j E + A W T jW + A S T j S + A N T jN + S lj 

(£>.45) 

where 

A Pj = ^ A nb + IttAxpAijp 

(£>.46) 



^ - C P (Re p / E c )Ma.x( — F* , 0 ) ; A* w = C p {R e ;/E c )Max(Fl, 0 ) 
A l N = C p (Re;/E e )Max(-FZ, 0 ) ; A* s = C p {Re* p / E c )U*x{Fl, 0 ) 

= [~(7tu uj + jtvVj + 7 tpPj + Jthhj) 

+ hBj(T B o ~ To) + hjj(T JO - f 0 )]A®pAy p 
+ {pfhu + h±)( Pje - Pjw )Ay p + 0Thv(p jn - p js ) Ax p 
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(DA7) 


(D. 48) 



